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Using a Kac-Moody current algebra with C/(l/l) x [/(l/l) graded symmetry, we describe a class 
of (possibly disordered) critical points in two spatial dimensions. The critical points are labelled 
by the triplets {l,m,kj), where I is an odd integer, m is an integer, and kj is real. For most 
such critical points, we show that there are infinite hierarchies of relevant operators with negative 
scaling dimensions. To interpret this result, we show that the line of critical points (1, 1, kj > 0) 
is realized by a field theory of massless Dirac fermions in the presence of U{N) vector gauge-like 
. . . static impurities. Along the disordered critical line {1,1, kj > 0), we find an infinite hierarchy of 

' relevant operators with negative scaling dimensions {Aq|g £ IN}, which are related to the disorder 

average over the g-th moment of the single-particle Green function. Those relevant operators can 
be induced by non-Gaussian moments of the probability distribution of a mass-like static disorder. 
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I. INTRODUCTION 



The plateau transition in the integer quantum Hall effect perhaps provides the most compelling experimental 
evidence in support of the scaling hypothesis in the theory of localization close to the mobility edge. Scaling with 
temperature, sample size, as well as dynamical scaling have been observed in two-dimensional degenerate electronic 



> 

■ systems subject to a strong magnetic field . (How strong the magnetic field must be is controlled by the ratio of the 
, magnetic length to the mean free path of the system.) Assuming that there exists one and only one relevant diverging 
' length scale ^ at the transition, the experimentally observed scaling is consistent with a power law divergence of f on 
' the energy separation from the center of the Landau band (the mobility edge) with scaling exponent I'cxp ~ i ■ "^^^ 
same exponent has also been measured in numerical simulations of the plateau transition for two-dimensional non- 
interacting electronic system in the presence of a static (quenched) weakly correlated disorder and a strong magnetic 
field [Q. It is tempting to interpret these experimental and numerical observations as evidences for the existence of 
a quantum critical point for a disordered non-interacting two-dimensional system at which, as a function of a strong 
magnetic field, a zero temperature continuous phase transition between two insulating quantum Hall states takes 
^ place. We take this point of view in this paper. 
Q ' The present understanding of the plateau transition from a theoretical point of view is quite unsatisfactory for two 
O basic reasons. First, it has not been possible to construct a field theory describing a second order phase transition 
^ [ governed by only one diverging length scale with scaling exponent = |. On the one hand, a semi-classical 
^ , picture based on the continuum model of percolation yields, in the infinite strong magnetic field limit, the critical 
' exponent i/^^j. = f @ • Although the classical continuum model of percolation does not reproduce the correct scaling 
exponent for the plateau transition, it inspired Chalker and Coddington to propose a network model incorporating in 
a simplified way quantum effects such as tunneling and quantum interferences ||^ . Numerical studies of the Chalker- 
Coddington (CK) network model suggest that the physics essential to the understanding of the plateau transition has 
been incorporated in it, since, among others, they predict the scaling exponent « | for the localization length 
in agreement with "first-principle" numerical simulations On the other hand, motivated by the mapping of the 
problem of the Anderson metal-insulator transition onto an effective non-linear a model, Pruisken and collaborators 
proposed a non-linear a model with the addition of a new topological term to describe criticality at the plateau 
transition Q . This approach has been partially successful in that it predicts the existence of critical points describing 
transitions between some insulating phases. However, it is not clear whether those critical points correspond to the 
plateau transition. Moreover, it has not been possible to calculate their critical properties. Further attempts to identify 
the effective field theory describing the critical properties of the plateau transition have consisted in approximating the 
replicated (supersymmetric) non-linear a model with topological term by a generalized replicated (supersymmetric) 
spin chain ||] ( Q), or in approximating the CK network model by a replicated Hubbard chain |^]. Unfortunately, 
all these effective field theories have remained quite intractable to analytical methods. The intrinsic difficulty in 
calculating the critical indices of the plateau transition is that this transition is between two insulating phases. This 
should be contrasted with the conventional Anderson metal-insulator transition for which critical properties can be 
calculated perturbatively from the metallic regime . All attempts to calculate the critical properties of an effective 
field theory presumed to describe the plateau transition are of a non-perturbative nature and have thus far failed. 
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Second, a highly non-trivial property of the critical field theory describing the plateau transition is that its operator 
content must be consistent with the observed phenomenon of multifractality [|lO|. Multifractality, in the context of 
the plateau transition, means that infinitely many universal exponents, which are not linearly related, are needed to 
characterize the scaling of the generalized inverse participation ratios at criticality. [| It reflects the fact that the 
inverse participation ratio is not a selfaveraging observable, whereas its logarithm is. Duplantier and Ludwig have 
shown how difficult it is to find from a stable Lagrangian field theory a critical point characterized by an infinite 
hierarchy of operators related to the phenomenon of multifractality [Tl| . There are some perturbative calculations in 
2 + e dimensions and asymptotically close to the mobility edge of disorder averaged generalized inverse participation 
ratios |p2[-p^. They are interpreted as evidences for multifractality jl^. It would be highly desirable to obtain 
non-perturbatively a hierarchy of scaling exponents for disorder averaged generalized inverse participation ratios, and 
to show how they are related to the universal scaling exponents observed numerically from the scaling with system 
size of the generalized inverse participation ratios for fixed realizations of the disorder. 

It was recently shown that a periodic potential can induce a continuous transition between two quantum Hall 
states or from a quantum Hall state to a Mott insulator in an anyonic system (at zero temperature) subject to a 
magnetic field |l6| , |l7| . It is then natural to probe the fate of such a (pure) quantum critical point upon perturbation 
by some weak disorder. In fact, such a procedure was used by Ludwig in a different context to show the existence of 
a new random fixed point for the two-dimensional g-state Potts model {q > 2) with weak quenched bond randomness 
psf . Ludwig also showed that the random fixed point for the two-dimensional g-state Potts model is characterized 
by an infinite hierarchy of relevant operators implying a complex scaling behavior related to multifractality A 
convenient model undergoing an integer quantum Hall effect in the absence of disorder was soon thereafter constructed 
by Ludwig et al. |20|. The critical properties of the pure system can be described by one species of massless Dirac 
fermions in 2 -I- 1 space-time. In the absence of disorder, the mass m of the Dirac fermions plays the role of a relevant 
coupling inducing, as it changes sign, a quantized jump in the transverse static conductivity. Ludwig at al. then 
studied the effect of three types of quenched (static) Gaussian randomness on the properties of the pure critical point: 
(?) an Abelian vector gauge-like impurity potential A^, /i = 1,2, (ii) an Abelian scalar-like impurity potential V, 
(iii) and a mass-like impurity potential M. The independent strengths of the three impurity potentials are measured 
by the variances y j^.j of the white noise probability distributions VIA], V[V], and 7^[M], respectively. To leading 
order in the impurities strength, the three randomness represent marginal perturbations to the pure critical point. In 
the presence of only one type of impurities, it was shown in Ref. pQ ] that, up to second order, g^j^y-j is marginally 
irrelevant (relevant), and that 5^ is exactly marginal to all orders, thus generating a line of critical points starting at 
the pure limit of the Dirac fermion model. The line > 0, gy = 0, gj^.j = of critical points is in fact multicritical 
since Gaussian randomness in V and M induce relevant perturbations [2C[| . Another important property along the 
multicritical line g^ > 0,gy = 0,gM = is that, after proper infrared and ultraviolet regularizations, the amplitude 
of the typical (with respect to the randomness) critical wave function in a large box (eigenstate obeying periodic 
boundary conditions with vanishing energy) is expected to be selfsimilar down to a microscopic length scale (the 
ultraviolet cutoff) and to scale like a multifractal measure with the system size (the infrared cutoff). Because the CK 
model can be viewed as a random tight-binding Hamiltonian, there must exist a lattice regularization of a field theory 
of random massless Dirac fermions describing the plateau transition. Ludwig et al. then conjectured that there is, 
in the three-dimensional coupling space {51^ vm > ^ stable fixed point which describes the critical properties of 
the plateau transition in the integer quantum Hall effect. However, if this fixed point exists, it cannot be reached 
perturbatively from any point of the multicritical line g^ > 0,gy = 0,(7^ = 0. 

In this paper, we study models of random massless Dirac fermions in two spatial dimensions. These effective 
field theories describe, in the presence of static disorder, the physics of the low energy and long wave-length single- 
particle excitations of quasi two-dimensional electronic systems which are characterized by isolated Fermi points. 
Examples for which models of random Dirac fermions are appropriate include degenerate semi-conductors |21 ] , two- 
dimensional graphite sheets , tight-binding Hamiltonians in the flux phase |Q , and dirty d-wave superconductors 
in two dimensions |24|j2^ ]. Models of random massless Dirac fermions are closely related to models on the effect 
of randomness on criticality in statistical mechanics (e.g., two-dimensional Ising |26| , ^7[ | and g-states random Potts 
model ]l8| , p^ ), since, at the energy of interest, the pure Dirac system is already critical. In the pure system and at the 
critical (vanishing) energy, the spatial decay of Green functions is algebraic. Our goal is threefold. First, the study 
of two-dimensional models of random Dirac fermions might yield a systematic classification of some new disordered 



^ This is to be contrasted with gap scaling in usual critical phenomena. An example of gap scaling occurs in the two- 
dimensional Ising model. If the scaling of the susceptibility for the two-dimensional Ising model is controlled by the exponent 
t;, then the scaling of the n-th power of the susceptibility is controlled by the exponent nrj. 
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critical points. Second, their study might shed hght on the field theoretical realization of complex scaling behavior 
(such as multifractality) not present in the pure limit. Third, it might be possible to infer some unconventional 
properties of the field theory describing the plateau transition from the properties of the field theory describing the 
multicritical points > 0, = 0, g^^j — 0, although, as we will see below, they must differ in an important way. 



A. Summary of results 



We will assume only one type of randomness: impurities interacting with N independent species of massless Dirac 
fermions like U{N) = C/(l) x SU{N) (static) vector gauge fields would do. The U{1) case was studied by Ludwig 
et al. in relation to the plateau transition [^. The SU{N) case was studied by Nersesyan et al. in relation to 
d-w&ve superconductivity |2^. Both groups used the replica trick to treat the disorder. We use instead the graded 
supersymmetric representation of single-particle Green functions to treat the U{1) x SU{N) vector gauge-like 
randomness (see also Ref. |^). The advantage of the graded supersymmetric representation is that with it the 
operator content at criticality is much easier to obtain than by using the replica trick. The disadvantage is that it 
is limited to non-interacting systems. In addition to the Abelian randomness € U{1) with Gaussian probability 
distribution of variance 5^, the non- Abelian randomness is represented by the static vector gauge field G SU{N) 
with Gaussian probability distribution of variance g^ . We make the important assumption that both and are 
topologically trivial. For example, in the Abelian case, this means that A^ is smooth everywhere and thus is not 
associated with a configuration of magnetic monopoles. 

We first prove that the generating function for the disorder average of any given product of m retarded and n 
advanced single-particle Green functions factorizes into m + n pieces which are each independent of disorder, provided 
all energy scales in the single-particle Green functions vanish, i.e., at the critical energy. This factorization is unique 
to vector gauge-like randomness. For example, it already fails for scalar-like V or mass-like M randomness. Because 
of this factorization, we can relate the operator content of the effective field theory for the averaged single-particle 
Green function at vanishing energy to the operator contents of averaged products of retarded and advanced single- 
particle Green functions at vanishing energies. Moreover, we show that, when all energy scales vanish, averaged 
products of single-particle Green functions are algebraic functions for arbitrary strength > of the Abelian 
randomness, and in the limits g^ — 0, 00 for the strength of the non- Abelian randomness. In other words, the 
generating functions for averaged products of single-particle Green functions describe field theories at criticality along 
the lines > 0,(7^ = 0,oo, whenever all energy scales vanish. It is the unconventional properties of these critical 
effective field theories that are at the heart of this paper. 

For the disordered critical points > 0, = 0, 00, we calculate the averaged local density of states asymptotically 
close to the critical energy. We reproduce the results of Ludwig et al. Q for the Abelian case, and Nersesyan et 
al. for the non- Abelian case It is noteworthy to notice that in the Abelian case, the averaged local density of 
states can be constant or even diverge for strong enough randomness (compare with and Since the order 

parameter for the averaged local density of states breaks a continuous symmetry of the effective field theory (the 
chiral symmetry), any finite or diverging averaged local density of states at criticality (vanishing energy) signifies the 
breakdown of the Mermin- Wagner theorem. The physical origins of the breakdown of the Mermin- Wagner theorem 
are two-fold. On the one hand, the pure system of massless Dirac fermions is critical. As a result there exists a 
continuous chiral symmetry in the effective field theory for averaged product of single-particle Green functions at 
vanishing energies [ p3pl| ]. On the other hand, the effective field theory for the disordered system is non- unitary, 
whereas the Mermin- Wagner theorem is proven for unitary field theories with compact symmetries. A power law 
dependency on energy for the averaged local density of states of disordered Dirac fermions is unusual compared to the 
generic case of a non-critical averaged single-particle Green functions in the problem of localization at the mobility 
edge such as, say, the plateau transition [ p2p^ ]. Note that there are other examples of diverging averaged local density 
of states in two-dimensional disordered systems [^,0 . 

The new result of this work is that we can calculate exactly the dominant scaling dimensions of arbitrary integer 
powers q of the local order parameter corresponding to the averaged local single-particle Green function. They are 
given by 

\ = <1 - - ^ i^' + N'l) ■ (1-1) 

This result is non-perturbative. It yields a spectrum of exponents which are not linearly related, form a concave 
set, are always relevant, and can be negative: w —q^ for large q. The physical meaning of Eq. (^]^) is that the 
local single-particle Green function, if considered as a random variable, is not selfaveraging, rather it is log-normal 



distributed at criticality. Eq. (1.1) was already obtained in Ref. [p4[ for the Abelian case. 
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The existence of an infinite hierarchy of operators which are relevant at criticahty calls into question the stability 
of the disordered critical points > 0,93 = 0,oo. As was pointed out by Ludwig et al. impurities inducing 
a scalar-like randomness ^ or a mass-like randomness M with Gaussian probability distributions generate relevant 
perturbations with scaling exponents Aj. They then conclude that the points > 0, — 0,oo are multicritical [^ . 
Our result, however, shows that the points g^ > 0, g^ — 0,oo (excluding the pure limit) are unstable in a very special 

way, namely that there are infinitely many relevant perturbations {^ITmIt^ G IN} which are generated by all non- 
Gaussian cumulants of the scalar- like or mass-like probability distributions. (For disorder with Gaussian probability 

distribution, only the variance, say g\ ^ y = g^ g v m-i non- vanishing.) Unless there exists a symmetry of the 
underlying lattice regularization of the random Dirac fermions which forbids randomness in the scalar-like or mass- 
like potentials, the search for the disordered critical points of random Dirac fermions must take place in an infinite 
dimensional disorder induced coupling space. 

It is instructive to contrast our results on the operator content for the disordered critical points > 0, 5^ = 0, oo to 
the 0{n) non-linear a model in 2+e dimensions for which there exists a non-trivial critical point {t* , h*) = {t^{e) > 0, 0) 
in a two-dimensional coupling space. Both the temperature t and the magnetic field h are relevant perturbations, 
and both have to be tuned simultaneously to reach the critical point. The fact that the critical point is unstable with 
respect to t does not mean that there is no physics associated with the h direction: correlation functions such as the 
magnetization near the critical point will depend on both t and h. In our model of random Dirac fermions, there are 
infinitely many relevant directions, each additional direction corresponding to a new class of randomness (characterized 
by the number of non- vanishing cumulants). If one were to compute some correlation function (numerically, for 

example) close enough to the, say, multicritical line > 0,g\/\,j = 0, the correlation function would depend on the 
form of the probability distribution for the disorder (all gyjj^.j) not just on its strength (gyjj^j)- The multicritical 
line of Ludwig et al. pO| ] is unstable in the sense that it is not possible to find in its vicinity a simple scaling form 
(i.e., depending on a finite number of universal exponents) for correlation functions. It is interesting to note that 
the non-trivial critical point [t* ,h*) of the 0{n) non- linear a model in 2 -I- e dimensions also has infinitely many 
relevant operators However, they break the 0{n) symmetry and cannot be realized physically by switching on 
an external magnetic field. Hence, they do not destabilize the critical poi nt (t*, h*). Final ly, t he validity of the 2 + e 
expansion itself has been questioned, due to delicate issues of convergence g^,^. In Eq. (1.1), the scaling exponents 



{Ag|q G IN} are obtained non-perturbatively, a highly non-trivial result if one attempts to derive it using the replica 
trick. 

We recall that in the theory of localization, the multifractal analysis consisting in first generating, for a given 
realization of the disorder, so-called box probabilities from the squared amplitude of a normalized energy eigenstate, 
and then performing a finite size scaling analysis on their q-th moment, yields a hierarchy of scaling exponents 
{T{q) — D(q)(q~' l)\q G IN} When the t((3')'s are not linearly related, the wave function is said to be multifractal. 
Multifractal wave functions are observed for energies close to the mobility edge in two spatial dimensions. More 
importantly, the T(g)'s are observed to be independent of the realization of disorder pO| . It is then tempting to 
conjecture that {i) the scaling exponents {T{q) = D{q){q — l)\q e IN} are universal, and [ii) can be calculated from 
some critical effective field theory describing the critical wave function in the thermodynamic limit. It turns out that 
we can relate the spectrum of scaling exponents {A^|g G IN} for powers of the local order parameter corresponding 
to the averaged local density of states to scaling of the averaged generalized inverse participation ratios introduced by 
Wegner This is so because the critical behavior of averaged products of advanced and retarded single-particle 
Green functions is controlled by the averaged single-particle Green function. Again random Dirac fermions differ in 
this respect from the conventional critical behavior at the mobility edge in the theory of localization, since in the 
latter it is the averaged product of one retarded and one advanced single-particle Green function, in other words the 
amplitude and not the phase of the single-particle Green function, which controls criticahty. 

We find that, along the critical line > 0,(/^ = 00, the scaling of the averaged generalized inverse participation 
ratios with respect to the energy separation from the critical energy and asymptotically close to it is controlled by 
the scaling exponents 

. ^.,„ . 2 - (£^ , ^) , e N. (1.2) 



where z is the dynamical exponent: 



TT 



The non-linear contribution in q to the right-hand-side of Eq. (1.2) comes from A^. Consequently, the function 



T*{q) = zw{q), when continued to real g, is non-linear, concave and unbounded from below for large q. The strict 
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concaveness of r*(g) reflects the non-selfaveraging nature of the squared amphtude of wave functions with energy 
close to the critical energy. The unboundness of r* (q) arises because wave functions have not been normalized. Rare 
events, i.e., realizations of the disorder such that amplitudes of the wave functions show "giant" spikes, can then 
dominate the averaging procedure for large q. 

We would like to stress that the function T*{q) only equals T{q) for q — and q = I. It is thus too naive to use the 
sole information contained in the hierarchy of local operators with scaling dimensions {A^jg e IN} (see Eq. (LI)) to 
calculate T{q) from the critical field theory. Moreover, it is far from obvious that the strict concavity of T*{q) implies 
the multifractal scaling of the probability measure constructed from the critical wave function for a given realization 
of the disorder, i.e., the strict concavity of T(g), as is implicitly assumed in Refs. [ p^ , p^pl|j20| ] . The important issue 
of how to calculate the universal scaling exponents {T{q)\q € N} within our critical field theory remains open, since 
the scaling exponents {A^jg e IN} arc not simply related to the T(g)'s. An attempt to relate the A^'s to the T(g)'s 
will be published elsewhere |38| . 

On the other hand, scaling with respect to spatial separation and effective system size of the averaged two-point 
correlation between the q-th generalized inverse participation ratio and the r-th generalized inverse participation ratio 
agrees in form with the one proposed for non-normalizable box-observables obeying multifractal scaling ||39| , [ll| , |l0| . 
Since the existence of negative dimensional operators is crucial to the former scaling, it is tempting to conjecture that 
multifractality of the critical wave function might be related to the existence of negative dimensional operators in the 
effective critical field theory. 

To uncover the nature of the critical field theories with 5^ > 0, 5^ = 0, c», we study their current algebras. Without 
loss of generality, we restrict ourselves to Dirac fermions with Abelian vector gauge-like disorder. We show that any 
disordered critical point i?^ > realizes a Virasoro algebra with vanishing central charge. Moreover, we show that 
this Virasoro algebra is obtained from currents obeying aJ7(l/l)xJ7(l/l) graded Kac-Moody algebra. The vanishing 
of the Virasoro central charge follows from using an effective partition function which is chosen to be unity in order to 
perform averages over the disorder. The Kac-Moody currents must obey a f/(l/l) x [/(!/!) graded algebra, since, on 
the one hand, the Hamiltonian for our random Dirac fermions can be interpreted as a Lagrangian in two-dimensional 
Euclidean space with U{1) x C/(l) chiral symmetry, and, on the other hand, bosonic copies of the Dirac fermions have 
been introduced to perform averages over disorder. After having introduced bosonic copies of Dirac fermions, the 
effective field theory cannot be unitary anymore. All unusual properties of the disordered critical point > can 
be traced to the bosonic sector of the effective field theory. 

Having identified the underlying current algebras for the disordered critical points > 0, we find new critical 
points by considering more general current algebras. We simply require that the critical points be described by 
two-dimensional Conformal Field Theories (CFT's) |^ with vanishing central charges and obeying C/(l/l) x t/(l/l) 
Kac-Moody current algebras. All such critical points are classified in terms of three numbers {l,m,kj). For odd 
integer I, the spin of Dirac spinors, which is ^ in the absence of disorder, remains half-integer as the U{1/1) x [/(l/l) 
symmetric interaction is switched on. For any real m and fc^ , the Virasoro algebra constructed from the currents 
generating the U{1/1) x [/(l/l) Kac-Moody algebra has vanishing central charge: a necessary condition to describe 
disordered criticality. Finally, the quantization of m follows from the requirement that the CFT be local. The Dirac 
fermions with Abelian vector gauge-like disorder provide a Lagrangian realization of the CFT with 



(/,m,fc^.) = (1,1,^) 



(1.4) 



We also construct Lagrangian realizations for all other CFT's. However, whether or not the remaining critical points 
describe criticality in systems of disordered non-interacting electrons is still an open question. 

For each critical point (Z, m, fc^), we obtain the operator content from the current algebra. We show, for example, 
that there exists an infinite hierarchy of local operators transforming like C/(l/l) x [/(l/l) singlets with the spectrum 
of scaling exponents 



(l,m,k-) 



mq + q^{l — 



g e IN. 



(1.5) 



Note that Eq. (1.1) is recovered for (Z, m, k^) given by Eq. (1.4). Conditions on /, to, and to insure that no operator 
with negative scaling dimensions be present at criticality are easy to obtain. When operators with negative scaling 
dimensions are present in the operator content and are not forbidden by some additional symmetry requirement, the 
critical point is unstable in the special way described above. 
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B. Outline 



The paper is organized as follows. We briefly review the supersymmetric representation of single-particle Green 
functions in the theory of localization in section M. The operator content of Dirac fermions in the presence of Abelian 



and non- Abelian vector gauge-like randomness is derived in sections III and IV , respectively. The scaling with energy 



of the averaged local density of states, averaged generalized inverse participation ratios, and of the averaged spatial 



correlations of generalized inverse participation ratios is derived in section |V|. In section VI, we study the algebraic 
structure of Conformal Field Theories related to Dirac fermions in the presence of vector gauge-like disorder. This 
section can be read independently from the previous ones. We conclude with a discussion of implications from our 
work for the plateau transition and the random XY model. In appendix ^ we give a detailed proof of the Sugawara 



V] 



construction of section |V|. Our last appendix ^ is devoted to Lagrangian realizations of the CFT's constructed in 
section 



II. AVERAGED PRODUCT OF SINGLE-PARTICLE GREEN FUNCTIONS 

We are concerned with calculating the disorder average over arbitrary products of matrix elements of the resolvent 

G{x,y,u;±w,V) = {x\ „ \, , . Iv), V > 0- (2.1) 

oj — Hq — V ±11] 

Here, Hq is any single-particle Hamiltonian in the absence of disorder, whereas V is any static single-particle potential 
representing disorder. Both are assumed to be local in space. The limit ^ 0+ gives the retarded and advanced 
(— ) single-particle Green functions in the energy representation for a given realization V of the disorder in a statistical 
ensemble with probability distribution VlV] . The coordinates x,y always denote spatial coordinates. They also 
implicitly denote degenerate internal degrees of freedom, unless explicitly stated. 
Impurity average over products of single-particle Green functions is defined by 

)({x„y„c^,±i,7}) ^ 



(2.2) 



\^^G{x,,y,,u;, +iri;V)j ^Y[G{x^^j,y^^j,uj^^^ - i7j;V)j = 

I V[V]V [V] [\{ G{x, ,y,,u,+ iTy; F ) j [\{ G(x„+,. , y„,+^ , c.,„+^. - ir,; V) 

We will be mostly interested in the averaged density of states per volume and per energy 

p{ijj) = lim -^7 tr [Gfx, x, w — iry; y) — G(x, w -f i?7; V^)], (2.3) 
))^o+ 27ri 

the impurity average over the q-th power of the smeared density of states per volume and per energy 



A(9'(w,ry) = {2T:i)-'i {tT:[G{x,x,uj-ir]-,V)-G{x,x,Lu±i7^]V)]}\ (2.4) 

and the impurity average over the product 

V'^^^'^\x,y,uj + i-q;y,x,uj - IT]) = G{x,y,uj + iT];V)G{y,x,uj-iiT,V). (2.5) 

The trace is performed over any internal degrees of freedom such as spin, orbital degeneracy, etc. Averaging over 
randomness is assumed to restore any spatial homogeneity which can be broken by a particular member V of the 
statistical ensemble. Hence, p{uj) and A^'^^uj,!]) are assumed not to depend on x. There exists an important Ward 
identity 



/ 



d'^x lim v\tTGix,y,uj + i7];V)\^ = T^Piod^), (2.6a) 



lim r] / d'^x \tr G{x,y,uj + iri;VW = irpiuj). (2.6b) 
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In Eq. (2.6a), the limit ry ^ 0+ is tak en be fore integrating over space. States which diffuse at infinity do not 
contribute to the right-han d-sid e of Eq. (^.6a ) for any finite separation \x — y\. These are the extended states. The 
only contributions in E q. (2.6a) come from localized states. On the other hand, both localized and extended states 
contribute in Eq. ( 2.6b|). W e will always be working in the thermodynamic limit, i.e., we will be concerned with the 
order of limits of Eq. (2.6b). 



A. Graded supersymmetric treatment of the disorder 



In this subsection, we recall how the representation of the single-particle Green functions in terms of a path integral 
over anticommuting (Grassmann) and commuting variables allows to perform the impurity average over the resolvent 
before taking its matrix elements |p8| . To this end, we simply use the Gaussian integrals 

dil;*dip e"'^*"'^ = a, Va complex, (2.7a) 

V'V'* e-'l'''"'' = 1, V6 complex, (2.7b) 

for the pair "0*1 "0 of Grassmann variables, and 

d(p*dip e^'^'^-'P ^ Re a > 0, (2.8a) 

a 

dip*d(p ip*ip e-'^''"^ " a / e"'^*"'^. Re a > 0, (2.8b) 

L J -^'a^ = l^e«>0' 2.8c 

J dtp* dtp e ^ a 

for the pair Lp* ,p> of complex variables and generalize them to Gaussian integrals for fcrmionic coherent states and 
bosonic coherent states, respectively. 

In this way, the denominator in the fermionic path integral representation 

can be rewritten as a Gaussian integrals over bosonic coherent states: 
G{x,y,u} ±\r];V) = 

i (v^Vi^ /py,*!?^ V(a;)V'*(2/) e+'/'''*'^*^*""^«"^^+"'l'''+'/'^*''*^''"^^«~^^+'''l'^. (2.9b) 



We could have equally well chosen the bosonic path integral representation 

G{x,y,uj±iri;V) = Ti^^ ^. r r,, .... , , (2.10a) 

which can be rewritten as 

G{x,y,uj±iiT,V) = 

Ti [v^*V^P (^(x)¥>*(y) e+'/^'[^^)^""''""''^+'''J'^+'/'"'[(^^(""''"-^)+"'l'^. (2.10b) 



Eqs. ( pJb| ) and ( pTOb| ) can be combined into 
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X e 



2 

+i / r[(±)(u~Ha-V)+ir,\-^+i J v*l{±){u;~Hg-V)+in]v ^2 H) 



Eqs. ( ^.9b| ,2.10b,2.11) are graded supersymmetric representations of the matrix elements of the resolvent. They require 



Hq + y to be quadratic in the fermionic {ip) or bosonic {(p) coherent states and thus can only be used when the pure 
system is non-interacting. 

In the graded supersymmetric representation, average over the statistical ensemble of impurity potential can be 
performed before calculating the matrix elements of the resolvent, and by so doing, averaged matrix elements of the 
resolvent are obtained from Green functions of an effective interacting theory for fermionic and bosonic coherent states. 
The effective interaction resulting from performing first the impurity average depends on two distinct properties of 
the impurities. One property is the nature of the coupling between the impurities and the single-particle eigenstates 
of Hq. It is the same for all impurity potentials in the statistical ensemble. The second property is the probability 
distribution V[V] for the impurity potential. The scaling hypothesis in the theory of localization assumes that given 
one type of impurity potential, there can exist a critical energy w^, the so-called mobility edge. Upon approaching 
this critical energy, the scaling hypothesis assumes the existence of one relevant diverging length scale. Whereas the 
existence of the mobility edge is determined by the nature of the coupling between the impurities and the single- 
particle eigenstates of Hq, the microscopic properties of the probability distribution V[V] should be unobservable 
close to the mobility edge. In other words, properties of V[V] on length scales much smaller than the diverging 
relevant length scale should not be observable. It is then natural to consider white noise probability distributions, 
i.e., statistical ensembles such that 



Vaix)V,{y) - v'6,,Six-y), (2.12) 

where a, b collectively denote degrees of freedom such as spin, orbital degeneracy, and so on, to describe a statistical 
ensemble with short-range spatial correlations. The energy scale v then alone characterizes the impurity strength. 
In the spirit of the scaling hypothesis, it is usually assumed that close to the mobility edge the choice of a Gaussian 
probability distribution, 

V[V] cx [|e-i^*^[^'(")], (2.13) 

can be made without loss of generality. Averaging over the impurity potentials then amounts to simple Gaussian 
integrals inducing an effective interaction which is quartic in the coherent fermionic and bosonic fields. In the graded 



supersymmetric representation, the average defined in Eq. (2.2) becomes the 2(m + n)-point Green function 



r(™'")({2:fe, y„ o^fe ± iT?}) = i"-" ^ Ijl {xMt{y^)j \ \{ An+,{Xm+,W.n+j{y,n+,)j j ■ (2.14a) 
Here, the expectation value is taken with respect to the partition function Z^^'""^ defined by the weight 



e'^o«'"'({-.±i'<}) ^ y e'^.""'+'«^"'-^r^"\ (2.14b) 

Sim) ^ + ^ / rfrf^ -H,-V) V^, + ^* (c., -H,-V) (^J , (2.14c) 

Z^l 

n „ 

Si"^ = - E / - ^0 - + ir , ^) ^ (^* , ^)] , (2.14d) 

Sim+n) ^ d'x [rk^ + ^l^,], (2.14e) 



fc=i 



(2.14f) 



and the measure 



V[r,i^,^*,^] = YlYl di;lix)dij,ix)dipUx)dif,{x). (2.14g) 



X fc = l 
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Notice that by construction -Z^g = 1- 

The action sij"^~^^^ insures the convergence of the path integral. It also has the largest symmetry, being diagonal 
in all degrees of freedom. In particular, it possesses a local U{m + n) symmetry where m and n represent the number 
of retarded and advanced single-particle Green functions over which the impurity average is performed. The action 
for the retarded sector has a smaller global U{m) symmetry when all energies Wj, i — I, ■ ■ ■ ,m are equal. The 

(n) 

same holds in the advanced sector S'a with n replacing to. Hence, in the special case 

^1 = • • • = = ^ri = • • • = ^m+n = ^a; (2-15) 

the effective action S'^^'"'' is symmetric under global U{m) x U{n) transformations acting on the index fc = 1, • • • , m+n 

of the coherent states. Another important property of the effective action S'^g '"^ is that when all energies are set to 
zero in the single-particle Green functions: 

= t^i = ••• = t^„+„, (2.16) 

then Sr"^^ + is invariant under the group U (to, n) which leaves 

diag ( +1,- - -, +1 , -1,- ■ ■,-! ) (2.17) 

m n 

unchanged. 



III. ABELIAN VECTOR GAUGE RANDOMNESS 

From now on, we restrict the generality of our discussion to two spatial dimensions and to the single-particle 
Hamiltonian describing one species of massless Dirac fermions in the presence of static Abelian vector gauge-like 
randomness pO[ |. In other words, the Hamiltonian of the pure system is the two by two matrix 

Ho = --^lA- (3.1) 

Here, /i = 1, 2 denote the two spatial directions and 7^^, 72 are any two of the three Pauli matrices, 75 = — i7i72 being 
the third. The coupling of the impurities to the single-particle states is through the dimcnsionfuU vector gauge-like 
potential A^{x): 

V{x) - -V9^A^{xh^. (3.2) 
The probability distribution for the random vector gauge-like disorder is Gaussian: 

r[A^,] cx e-*/'^'^^'("). (3.3) 

In the graded supersymmetric path integral representation of the Green functions, the fcrmionic coherent states are 
denoted (ipji/j), while the bosonic coherent states are denoted {(p^^p)- Both fcrmionic and bosonic coherent state s 
transform like Dirac spinors in two-dimensional Euclidean space, and thus implicitly carry spinor indices. Eq. (2.14) 
becomes 



r("^")({a:„y,,^, ±iry}) = i"-'" ^ \X{i^AxMy^)j j^n^™+.(^™+.)^"+.(2/m+,)j y > (3.4a) 

= + E / {^^ - i^^A^) + Lo,] + V) ^ (^, V)} , (3.4c) 
1=1 •' 

n „ 

^i"^ = - E / [i^A'^^M - iV^I^^) + ^rn+j + i^, V-) ^ (<^, <P)} , (3.4d) 

Sim+n) ^ rjJ2 (^.V-fc + ^kfk) , (3.4e) 
k=i •' 

m+n 



fc=l 
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Integrating over the random vector gauge-like potential induces a quartic coupling mixing all the indices labelling 
the m + n single-particle Green functions. It would therefore appear that the effective action S'^^"*'"^ is hopelessly 
complicated. This is not so however in a very special case, namely when all energy scales are set to zero: 



m+n 



= 0, V = 0. 



(3.5) 



Indeed it is then possible to perform a very simple trick to decouple entirely the Dirac spinors from the random 
vector gauge-like potential. To see this, we first make use of the fact that, since x belongs to a two-dimensional 
manifold taken to be the Euclidean plane and since the vector potential {A-^^x) , A2{x)) has only two components, 
can always be decomposed into a transversal ^i{x) and a longitudinal ^2{^) component, provided is free of 
singularities associated with configurations of magnetic monopoles. Accordingly, we can use the Hodge decomposition 



e d 



(3.6) 



We have used the completely antisymmetric second rank (Levi-Cevita) tensor e^^ with — 1- Parametrizing the 
random vector gauge-like fields in terms of their transversal and longitudinal components induces two changes in their 
probability distribution. The first one comes from the Jacobian induced by the change of variables. The second one 
comes from the change in their probability distribution. We thus have |42 



PK] ^K] = Dct [d^] 



(3.7) 



We are implicitly assuming that the vector gauge fields vanish at infinity. On the last line, the Jacobian of the 
transformation is rewritten in terms of a path integral over fermionic coherent states {a, a). We will refer to these 
degrees of freedom as ghost fields. They do not play an important role for an Abelian vector gauge impurity potential, 
but they must be accounted for in the case of non- Abelian vector gauge impurity pot entials as we will see later on. 
The reason for which Abelian ghosts are not important is that the Jacobian in Eq. (3.7) is independent of the Abelian 
vector gauge-like configuration A^. 



After the change of variables Eq. (3.6), the minimal coupling between the Dirac spinors and the impurity potential 
takes the form 



in the fermionic coherent state sector, and 



fkln 



fc = 1,' 



, m + 71, 



(3.8a) 



(3.8b) 



in the bosonic coherent state sector. Again we make use of the fact that we are effectively working with covariant 
gauge fields defined on a two-dimensional manifold so that the important identity 

7^ 75 = 7^> = 1,2, (3.9) 

holds. It is then easy to verify that for the fc-th Green function in the product r^'"'"^ the change of variables 



'fk 



Vk e 



Vk 



"Pk- 



(3.10a) 
(3.10b) 



decouples the Abelian vector gauge-like fields from the Dirac spinors [l43| . This local transformation can be repeated 
for all TO -|- n single-particle Green functions. It is not a pure gauge transformation since it involves the longitudinal 
component <f>]^. It is the product of a pure gauge transformation (the exponent $2) with an axial transformation (the 

exponent ^i). 

The axial factor of the decoupling transformation Eq. ( [3.10| ) is non-trivial in many respects. First, in contrast to 
the pure gauge factor, it does not leave the terms proportional to the energy scales w^., 77 unchanged: 



1, • • • , TO -I- n. 



(3.11) 



Second, it does induce a non-trivial Jacobian for both the fermionic and bosonic coherent state measures, as long as a 
fully gau ge i nvariant regularization of these measures is used [ p4[ . However, the Jacobian induced from the fermionic 
measure fEsI 
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(3.12) 



cancels that from the bosonic measure (this cancellation does not hold for non-Abelian gauge fields due to the ghost 
sector). Finally, the manifold defined by all axial transformations is not compact. 

The impurity average of any product of single-particle Green functions factorizes when all their energies are set to 
zero: 

hm r(™-")({x„y„±i,7}) = 

?7— >0+ 



)'^m+jiym+j) 



e75*i(^>)e75*i(?/i) 



J^e+'*2(a;,)g-ii'2(a,) 



n 



(3.13a) 

(3.13b) 
(3.13c) 

(3.13d) 

(3.13e) 
(3.13f) 
(3.13g) 

Thus, the impurity average of any product of single-particle Green functions can be obtained from the correlation 
functions generated by the corresponding critical Euclidean effective action 



The expectation value with respect to the spinor coherent states is |^ 

m „ 

71 „ 

i=i 

The expectation value with respect to the longitudinal component of the Abelian gauge potential is 
The expectation value with respect to the the transversal component of the Abelian gauge potential is 
Finally, there is an innocuous average over the ghost coherent states 



(3.14) 



The impurity strength does not appear explicitly, being hidden in the rotation from the original fields to the new 
ones. The effective action is criti cal in the sense that all correlation functions are algebraic as we will see below. 
From the factorization in Eq. ( p. 13 ) follows immediately that given m -\- n pairs of points x^. , yj. on the Euclidean 

plane, the 2(TO-|-n)-point correlation function T^"^'^\{xj., y^,, 0}) can be calculated by using any effective action S'cf'''* 
as long as p > TO and q > n. For example, the average local density of states lim^^^o p(^) can be calculated from all 
effective actions S'if'*'' with cither p > 1 or g > 1. We now turn to the operator content of the effective actions S'c™'"''. 



O f cours e, integration over the bosonic coherent states is not guaranteed to converge anymore having set t] to zero. Hence, 
Eq. (3.13b) is only a formal expression unless it is itself regularized. 
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A. Operator content at criticality for Abelian vector gauge randomness 



To study the operator content of the critical action in Eq. (3.14), it is useful to introduce the so-called chiral basis 
for Dirac spinors V', V': The chiral basis is defined by 



(1±75) V', 



= (^7i, (^^t = <y="^^(l±75): 'P± = ^(1±75) V- 



In the chiral basis, the decoupling transformation, Eq. ( 3.10 ), becomes 



t 't 



i>k± 



Vk±, 



/fc = 1, 



1, 



, TO + n, 

, m + n. 



(3.15a) 
(3.15b) 



(3.16a) 
(3.16b) 



The terms in Eq. (3.11) which couple to the energy in denominators of single-particle Green functions are represented 
by 



'^k'fk 



^fc'+e"'^*^^L + V'fc-e+2v^*ii/;^+, A: = 1, • • • , to -f n. 



''I 



(3.17a) 
(3.17b) 



provided (71,72:75) sse chosen to be the usual three Pauli matrices. The effective actions in the retarded and 
advanced Dirac spinor sectors take a simple form in the chiral basis. To see this, we first think of the Euclidean plane 
{x ~ (xi, X2)\xi, X2 G IR.} as being a real section of the two-dimensional complex manifold with local coordinates 



^2: 



(3.18) 



where Xi and X2 are now arbitrary complex numbers. The effective actions in Eqs. (3.13c, 3.13d) now become 

m 

= +2i^ d, + i^'l 9, + d, + 9, ^:_) , (3.19a) 

n 

i=i 

respectively. The advantage of the chiral basis is that the two-point correlation functions in the decoupled spinor 
sectors are very simple. The only non-vanishing ones in the retarded sector are (fc, Z = 1, • • • , m): 



(^^+(z,zX(0,0)) ,^ , = (^^^(z,z)^;t40,0)) . , = ^' 

In the advanced sector, they are the same except for a relative sign [k,l = to -I- 1, • • • , to -I- n): 

(v''fe+(z,z) v^;Uo,o)) ,^ , =(^',+(z,z)^;;(o,o)) ,^ , 



(^l._(z,z) V;l(0,0)) ,^ , =(^'fe_(z,z)^;i(0,0)) , 



^kl 
2-KZ' 



(3.20a) 
(3.20b) 

(3.21a) 
(3.21b) 



Two-point functions between retarded and advanc ed spinor fie lds vanish. In the chiral basis, t he two-point correlation 
functions are thus either holomorphic as in Eqs. (3.20a 3.21a) or antiholomorphic as in Eqs. (3.201 ,3.21b). 
Another important set of two-point correlation functions that we will need is 



(3.22a) 



where the conformal weights h and h are given by 



12 



Stt' 



(3.22b) 



The proportionality constant on the right-hand-side of Eq. ( 3.22a ) depends on the short distance cutoff used to 
regulate the free scalar field theory. The parameter / is analytically continued from the real numbers to arbitrary 
complex numbers. In particular, we will be interested in purely imaginary values of / owing to the non-compact 
nature of the decoupling transformation in the longitudianal (^i) vector gauge sector. 

For any given Euclidean effective action Sc™'"' in Eq. ( 3.14| ), we can calculate the scaling dimensions of all possible 
local operators, thus explicitly showing that all correlation functions are algebraic. To illustrate how this is done, we 
begin with the effective action S'ct'"^ for the averaged retarded single-particle Green function. One verifies that 



>+(z,z)VT_(0,0) 
;^_(z,z)^L(0,0) 
as follows immediately from, say, 

;vI(^,^)V'±(o,o) 



¥>_(z,z)^U0,0) 



1 

2^' 
1 

2^' 



(3.23a) 
(3.23b) 



>2(^,^)^±(o,o) 

X ^gi(±iV5l)*i(z,«)gi(Tiv^)*i(0,0) 



The bracket denotes an average performed with the effective action 



(3.24) 



(3.25) 



We see that any impurity (5^) dependency coming from averaging over precisely cancels that coming from averaging 
over $2- Hence the scaling properties of the fields ipj. and ipj. are not affected by impurity averaging. This is not so, 
however, when calculating the two-point function for the composite operator 



0(z,z) = ^Tp- ■0+ "0- V+ V- ) (-Zj -z); 

which is defined through point-splitting {rrii, • • • , mg are arbitrary positive integers). Indeed, 

X {^T^ (z, -z)^'!"^^ (0, 0) ) ( ^ (z, z)^r ^ (0, 0) ) 

X (^'i-^ (z, z)^'r« (0, 0)) ^ (z, z)^T' (0, 0)) 

X (v.'"^(z,z-)^'i"«(0,0)) (^r-^(z,z)^r^(0,0)) 



(3.26) 



(3.27) 



The conformal weight h (h) of ''/'+(_) scales quadratically with m in contrast to that of which scales linearly 

with m. This is perhaps most easily seen with the help of Abelian bosonization when calculating correlation fun ction s 
of V'i's. Exponentials in the longitudinal and transversal components of the random gauge fields scale as in Eq. ( 3.22| ) 
with 



/i = — i (rrii — ~ ™3 + '^4 + '^5 ^ '^6 ~ ™7 + ^s) ■ 

/2 = + (tox + "^2 ^ ^'^S ^ "^4 + ™5 + ™6 ~ ™7 ~ '^s) J 



(3.28a) 
(3.28b) 



respectively. Notice the difference in the relative sign of the m's in and /j. We thus conclude |34| that O scales 
like 
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2h--2h 
Z . 



where the conformal weights are given by 

h 
h 



1 - 

2 . 
1 r 
2 



rrin 



171. 



Too 



{fl 



l/l 
l/l 



|2\ 9a_ 
' 47r 

' 47r 



(3.29a) 

(3.29b) 
(3.29c) 



Irapurity averaging does not change the scahng dimensions of local operators when /| — \fi\'^. On the other hand, 
the non-compactness of the decoupling transformation allows for the possibility that local operators have negative 
scaling dimensions when \fi\'^ ^ /| . 

Negative dimensional operators are now easily obtained in the theory with effective action S'^^'^\ For example, 
consider the composite operator „ defined by 



n-, 



n-^ > 0, > 0, 

n-^ < 0, rij > 0, 

> 0, 77.2 < 0, 

n-^ < 0, < '^j 



(3.30a) 



and use Eqs. ( 3.2^ 3.29). The conformal weights of these operators are 



2tt 



2tt 



'-l"'2- 



(3.30b) 



This demonstrates that for any given value of g^, there are infinitely many local composite operators with negative 
conformal weights. This is very different from unitary two-dimensional conformal field theory, say the U{\) Thirring 
model |4^] (the action for the Euclidean U{1) Thirring model differs from S'ir'^'' only in that the ip sector is absent), 
which cannot support operators with negative dimensions. Operators like are not mere curiosities. For one 

thing, they are generated in the effective graded supersymmetric model if non-Gaussian moments in the probability 
distribution of a mass-like random perturbation of the Dirac Hamiltonian are present. Furthermore, they yield direct 
information on the impurity average of products of single-particle Green functions as we now show. 

The generalization to the calculation of the scaling dimension of any local operator in the theory with effective 
action S, 



{ra,n) 



is now obvious. The most general local composite operator 

o{z,-z) = n (<'^^v^l"^^c+^^C-^vi?vi"-^r+-^r-!^ 
fe=i 

where point-splitting is implied and all to's are positive integers, scales like 
(Ot(z,z) 0(0,0)) 

/l = - > iTTlTl- + TTlJj, -h TO 



{z,z), 



oc z 



-2h--2h 
Z 1 



h = 



fc=l 

2 H [^2fe - 
fe=l 

m-t-n 

k=l 
m-\-n 



5k 



'■6k 



'■Ski 



(/I -l/l I 
(/I -l/l I 



2^ 9a 
Stt' 

2^ 9a 
Stt' 



''2k 



m^y. + TO 



4fe 



'Sfe 



TO 



6fe 



TO. 



7k 



''2k 



'4fe 



''6k 



TO 



7k 



mak)- 



k=l 



An important consequence which follows from Eqs. dpl-llll) and Eq. (^isil) is that 

(Ot(z,z) 0(0, 0))^(i.o) ^ (7't(z,z) 7'(0,0))^(™^,.,) , j = 0,---,to. 



(3.31a) 

(3.31b) 
(3.31c) 

(3.31d) 

(3.31e) 

(3.31f) 

(3.32a) 
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where 



m~ m« m-. mo 



yj=l j = l k=l 1 = 1 



(3.32b) 

(3.32c) 
(3.32d) 



Eq. (3.32a ) relat es correlations in a model with several copies of the pair of spinors (ipj^p) to a model with a single 
pair. Eq. ( 3.32a ) is only strictly true in the limit z,z-^0. There are subleading correctio ns to this algebraic decay, 
in particular when j ^ 0,to. It is an example of an operator product expansion. Eq. (3.32) allows to relate the 
operator content in the theory with effective action Scl'^'^ to impurity averaging of higher powers of the single-particle 
Green functions. For example, it tells us that the scaling dimension of 5'_„j,„, m a positive integer, in the theory 

with effective action Scl'^^ is the same as the scaling dimension of the operator whose expectation value in the theory 
with effective action S'c™'"'' yields the impurity average [tr G{x, x, 0; y)]™. Indeed, we recall that 



[tr G{x,x,Q]V)Y 



\fe=l 



{x) 



(3.33) 



After expanding the product on t he rig ht-hand-side of Eq. (3.33), we can calculate the scaling dimensions of any of 
the resulting 4™ terms using Eq. (3.31). They all share the scaling dimension of 



m 



k=l 



which in turn can be extracted from the algebraic decay of 

{V\x)V (y))5(™,o, cx ((¥._^;)"(a;)(^+^i)"(y))^,,,„,. 



(3.34a) 



(3.34b) 



IV. NON-ABELIAN VECTOR GAUGE RANDOMNESS 



Tight-binding one-body Hamiltonians, for which the Fermi surface is made of a collection of isolated points in the 
Brillouin zone, are well described by, say, N species of Dirac fermions, if one is only interested in the long wavelength 
and low energy properties of the system. In general, impurities induce one-body interactions among the N different 
species. In this section, we consider such a two-dimensional system in which impurities can be represented by a static 
Abelian vector gauge-like interaction on the one hand, and by a static non-Abelian vector gauge-like interaction on 
the other hand. The Abelian vector gauge potential does not induce scattering events between different species. 
The non-Abelian vector gauge potential does. The index /i = 1, 2 denotes the two spatial components of the gauge 
fields, and a = 1, • • • , N"^ ~ 1 labels any basis {T°-} of traceless Hermitean N x N matrices. The Hamitonian of the 
pure system is the 2A^ x 2A^ matrix 

^^0 = -i^7^5^, (4.1) 
where / is the N x N unit matrix. The impurity potential is given by the 2A^ x 2A'^ matrix 

V{x) = 7^ + B^X^)] , (4.2) 

where 

= B; (4.3) 
transforms like the adjoint of the Lie algebra su{N). Our conventions for the generators of su{N) are 

tr(r"r'') = [T'',T''] = ir''"T", f^'^dfcd ^ ]^§ab^ ^jj^^ ^ l,...,iV2-l. (4.4) 
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Any impurity average over a product of single-particle Green functions can be obtained from an effective partition 
function. We are only interested in the single-particle Green functions in which all energies have been set to zero. 
Without loss of generality, we will concentrate on the averaged single-particle Green function. It can be obtained 
from the expectation value of 'ip{x)'>p{y) + tp{x)ip{y) with respect to the partition fimction 

Z ^ J / J ViB^] e'S"'^^^ (4.5a) 



A^A^ + — tr (B,.B„) ^ (4,5b) 



The Dirac spinor -0 denotes a Grassmann coherent state, a bosonic one. Both transform like the fundamental 
representation of su{N). The vector gauge fields are still dimensionfuU, although they have been rescaled compared 



to section [II. The positive impurity strengths and act like inverse masses for the Abelian and non-Abelian 



vector gauge fields, respectively. We show in this section that the partition function in Eq. (4.5a) factorizes into 
four independent sectors. However, in contrast to the case of Abelian vector gauge impurity only, we show that Z 
does not describe a theory at criticality for all but one value of > 0. Criticality is proved for the special case 



g^ — oo. In this limit, C in Eq. (4.5b) is shown to be conformally invariant, and to yield an energy-momentum tensor 
with a vanishing central charge in its operator product expansion in agreement with Ref. ]4^ . The vanishing of the 
central charge follows from the fact that by construction Z = 1, and thus serves as a consistency check. The operator 
content of Z is constructed for g^ = 0, g^ = oo and we obtain the new result that an infinite hierarchy of negative 
dimensional operators is present at the (non-Abelian) disordered critical point. 



A. Operator content at criticality for non-Abelian vector gauge randomness 

We begin by showing that it is possible to decouple the impurities from the spinor coherent states. To this end, let 

A, = A^-iA^, A, = A^+iAr^, (4.6a) 

B, = B^-iB^, B, = B^+iB^. (4.6b) 

The antiholomorphic components A^ and B^ of the gauge fields are independent from the holomorphic components 
A^ and B^ . Both antiholomorphic and holomorphic components of the gauge fields belong to the complex extensions 
u'^*{N) and u'^{N), respectively, of the real Lie Alebra u(N). The antiholomorphic and holomorphic components are 
a convenient choice when working with the chiral basis for the spinors since 

C = ijli2d, + iA,I + iB,)^j^ + ^pl{2d,+iAJ + iBj^p^ 
+ {2d, + iA-J + iB,) + <pl {2d, + lAJ + iSj 

-f ^A+A_ + —ir{B+B_). (4.7) 
"^9 A 9b 

We parametrize the components of the gauge fields in the Lie algebras u'^*{N) and u'^{N) by fields in the complex 
extensions U''*{N) and {/"{N) of the Lie group U{N): 

A-J = +i{2d,H,)Hi\ H, = e-'^*s = $2-i*i, (4.8a) 

AJ = +i{2d,HjH-\ = e-'^*s = $2+i*i, (4.8b) 

B, ^ +\{2d-fi,)G-,\ G,e SU"{N), (4.8c) 

B, = +i{2d,G,)G-\ G,e SU'{N). (4.8d) 

This change of variables costs the Jacobian 

Jgi, = Dot {2dg)I)et (2ajDet (V2)Det (VJ, (4.9) 

where Vj and are the covariant derivatives in the adjoint representations of SU''*{N) and SU''{N), respectively. 
For example 

SA, = i 5 {{28, G,) GJ') 

= {2d, . + i [i {2d, G,) Gji , . ] } (i {5 G,) GJ^) . (4.10) 
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J:)et{2dg) = 




,4 




d^x 2tr(/3!J_T'' 


i2ajT°a° 


Det(2aj = 








d^a; 2tr(/3lT'' 


i2a^T''a5i 


Det(V2) = 








rf^K 2tr(/3^T" 


iVi;T''a';_) 


Det(V,) = 


1 1? 


«-] 






iV^T''a^) 



We represent this Jacobian as a path integral over auxihary Grassmann fields (ghosts): 

(4.11a) 
(4.11b) 
(4.11c) 
(4.11d) 

The unit matrix / has been rescaled to T'^ = } I for notational convenience. 

We are now ready to decouple the vector gauge fields from all other fields. For example, the chiral transformation 

= iP'l G^'Hi\ V+ = H,G, (4.12a) 
V'l - ^'1 G-'H;\ ij_ = H,G, C (4.12b) 

decouples the V' sector from the vector gauge fields at the price of the Jacobian Q] 

_ Pet [i7^ {d^ + iA^j)] Pet [17^ {d^ + iB^)] 

Pet(i7,5j Det(i7,a,) ' ^^'''^ 

Similarly, the chiral transformation 

parj^a ^ Q_ 0arj.a q-1 ^ ^aya ^ q_ ^,aj,a q-1 ^ (4.14a) 

PIT" = G, P'^T" G^\ atT" = G, a'^T" G-\ (4.14b) 
decouples the ghosts from the vector gauge fields at the price of the Jacobian 

The e xpo nent 2N is twice the quadratic Casimir eigenvalue = iV in the adjoint representation as is apparent from 
Eq. (4.4). Finally, by performing exactly the same chiral transformation on the ip^s as we did on the ^/''s, we can 
entirely decouple the spinor sector from the vector gauge sector. Pecoupling the (^'s from the gauge fields costs a 
Jacobian 

= (4.16) 

which cancels the one coming from the sector. We are then left with the Jacobian jT^i ... from the ghost sector. 

Polyakov and Wiegmann have calculated jT^i ^jv. It is given by the exponential of the Euclidean Wess-Zumino- 
Witten action [Q: 

= e"(-2^)i^[«-~^^^-l, (4.17) 
where the Euclidean Wess-Zumino-Wittcn action is given by 

r[G] = ^ / dx.dx^ tr [{d^ G ) {d^ G-')] 

- / dxidx^dxg 6^,, tr [{d^ G)G-' [d, G)G-' {8^ G)G-' ] . (4.18) 



Before collecting all contributions to the effective Lagrangian density, we rewrite the probability distributions for 
the vector gauge fields in a more useful form. In the Abelian sector, the Gaussian weight is 

+ = "2^^^ [{2d,H,)H^'{2d,HjH-']. (4.19) 
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In the non-Abelian sector, the Gaussian weight is 



9b 9b 
With the help of the Polyakov-Wiegmann identity | 



^ tT{B,Bj = tr [(2a,G,)Gri(2a,GJG;i] 



r[G-^H] = T[G-^]+T[H] + ^ J dx^dx^ tr [{2d^G )G-^ {2d-^H )H-^] , 
the Gaussian weights for the Abehan and non-AbeHan vector gauge fields are rewritten 



1 

25^ 



9aN 



and 



-Itr {B,B,) = {r[G-'G,] - T[G-'] - r[G,]) 

9b 9b 



(4.20) 



(4.21) 



(4.22) 



(4.23) 



respectively. In the Abelian case, the mass term separates into two independent sectors. In the non-Abelian case, due 
to the topological term in the Wess-Zumino-Witten action, the mass term does not separate into two independent 
sectors. 

We conclude that the trace of the averaged single-particle Green function in which all energies have been set to 
zero can be obtained from the expectation value of 

tr [ilj{x)tp{y) + Lp{x)'f{y)] = 
tr 

tr 



(i7,G,^;) (x) [^/iG-'H-') (y) + {H,Gjy_) (x) [^'iGi'Hr') (y) 
(i/,G,^V) (x) [^'iG-'H-') (y) + {H,G^^'_) {x) [^'Ig^'H^') (y) 



with respect to the partition function 



Z 

S, = 
S, = 



N 

2 ^ d?x tr {l3'^T''2d-^a'lT^ + l3'^T''2d,a'^T'') , 

^ ) n ^ 



Ss = 



S4 = 



a,b=0 



9aN 



nHzH.h 



47r 

h 27V 

9b 



1 47r 1 47r 

r[G;iG,] + — r[G;i] + -r[G,]. 

9b 9b 



(4.24) 

(4.25a) 
(4.25b) 

(4.25c) 
(4.25d) 
(4.25e) 



We can use the invariance of the Haar measure of U''*{N) {W^iN)) under left and right multiplication to reparametrize 
the Abelian vector gauge components Hg and as H^^H^ and H^Hg. The partition function has factorized into 
Jour independent sectors. The same factorization takes place for partition functions representing the impurity average 
over an arbitrary product of single-particle Green functions, as long as all energies in the Green functions are set to 
zero. However, the factorization property does not guaranty criticality, i.e., algebraic decay of all correlation functions. 

The first sector describes a free massless relativistic theory with U{N/N) graded symmetry. All Dirac spinors 
transform like the fundamental representation of U{N). This theory is conformally invariant with primary fields 
and ipj. carrying the conformal weights 



ih,h) = (^,0) or (0,i), 



(4.26) 



depending on their chirality. The first sector has vanishing central charge 
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ci =0, 



(4.27) 



since the partition function restricted to this sector is unity. 

The second sector describes free massless relativistic ghosts transforming hke the adjoint of U{N). This theory 
has a conformally invariant energy-momentum tensor. The operator product expansion of the energy-momentum 
components yields the negative central charge p8| 

C2 = -2dim U{N) = -2N'^. (4.28) 

The third sector describes two independent and free scalar field theories and is thus conformally invariant. One has 
negative definite kinetic energy. Consequently, its primary field H^^Hg oc e'(^'*i^^ carries negative conformal weights 

h ^ h = (4.29) 

OTT 

The scalar field theory with negative definite kinetic energy corresponds to the longitudinal component of the Abelian 
vector gauge field. The second scalar field theory has the primary field H^Hg (x e^'*^^ with positive conformal weight 

h = h = (4.30) 

It corresponds to the transversal component of the Abelian vector gauge field. The central charge in the third sector 
is 

C3 = 2. (4.31) 

The last sector describes the non-Abelian vector gauge fields. In contrast to the Abelian vector gauge sector, it is 
not conformally invariant for all values of the non-Abelian impurity strangth g^ . Only the limit ^ oo ^ yields 
a conformally invariant theory, namely the SU{N)f. Wess-Zumino-Witten theory where k = ~2N is the level of the 
underlying Kac-Moody algebra j5l]. The primary field Gj^Gg carries the conformal weight ^l[] 

1 N'^ -1 

h ^ h ^ — T——, k = -2N, (4.32) 

N + k 2N ' ^ ' 

which is thus negative. The associated central charge [|l| 

is positive. In the limit — > oo , there exists a local non-Abelian gauge symmetry with G^Gg a pure gauge 
Integration over the pure gauge enforces the local constraint 

= 0, M=l,2, a=l,---,N^-l, (4.34) 

for the non-Abelian currents 

= i'T'^lf^i^ + ^T'^lf^V, /i = l,2, a=l,---,iV2-l. (4.35) 
Thus, only SU{N) gauge singlets can acquire an expectation value in this li mit. 



For all correlation functions calculated with the partition function in Eq. (4.5a) to be algebraic, one must take the 
limit gg — s- oo. In this limit, the existence of a local gauge symmetry requires a careful definition of the measure of 
the vector gauge fields. Gauge fixing is in fact implicit in the chiral transformation which decouples the non-Abelian 
vector gauge fields from the spinors and ghosts. Indeed, one should think of the chiral transformation as the product 
of a pure gauge transformation, which amounts to gauge fixing, and of an axial transformation to achieve complete 
decoupling |5^. After gauge fixing, only gauge inequivalent configurations must be summed over in the partition 
function. Hence, the partition function restricted to the non-Abelian vector gauge sector is given by 



^ It is amusing to note that for negative values (and thus unphysical) of the non-Abelian impurity strength gg = — the 
non-abelian sector is also "critical". 
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= J V[G'] e-(-2iv)r[G']^ (4 3g) 

in the limit oo. To verify the correctness of Eq. (4.36), we calculate the total central charge: 

c = Ci + C2 + Cg + Csy(A,)_^^ =0. (4.37) 

It vanishes as it should be. 

The conformal weights of all primary field {0;} of the SU{N)f. Wess-Zumino Witten theory have been calculated 
by Knizhnik and Zamoklodchikov . They are given by 

h = h = (4.38) 

c^ + k' c^ + k 

In the numerator, q (cj) stands for the quadratic Casimir eigenvalue of the irreducible representation / (l). Here, 
the representation / (/) is defined by the transformation law obeyed by 0; under SU (N) rotations generated by the 
holomorphic (antiholomorphic) component of the level k Kac-Moody currents. Thus, by analytical continuation of k 
to the value k — —2N, we see that all primary fields originating from the non-Abelian vector gauge impurity in the 
limit = oo have negative conformal weights 

The eigenvalues {c; } are known. Examples are 

_ (A^-2)(A^ + 1) _ ^ _ ^ ^ _ {N + 2){N-l) 

in the antisymmetric °, fundamental □, adjoint, and symmetric □□ representations, respectively. Here, □ represents 
a box in a Young tableau. We recall that all irreducible representations of U{N) and SU{N) can be labelled by Young 
tableaus made of N rows with boxes in the i-th row, provided > /2 > • • • > Jn-i > In — ^- For SU{N), 
fj^ always vanishes. It is known that the SU{N) irreducible representation with Yoimg tableau {/j } made of 

/ = ^iLi fi boxes, has the quadratic Casimir eigenvalue 

1 ^ f2 
= 2E[/' + (^ + l-20/.] - ^- (4.41) 
1=1 

We are now ready to compute the most relevant scaling dimensions controlling the expectation value of the local 
density operator ipip + ipip and of higher powers thereof in the limit g^ = 0, g^ —^ oo. Without loss of generality, we 

only need to consider powers of First and in agreement with Refs. |25|j47| ], the conformal weights of 

(ipl^^) (z,z) = (^'l G' ^V) (z,z), (4.42) 

are given by the quadratic Casimir operator in the fundamental representation □ of SU{N): 

'1 I N^-l\ 1 , , 

(4.43) 



2 N + k 2N J k=-2N 2iV2" 



Next, higher powers of (pL^fi^ are defined through point-splitting: 

rn m 

" {z, -z) = Hm n (^1^+) (z, , z- ) = Urn [] (^'!gVV) {z, , ^J- (4.44) 



We have already dealt with the product of Lp'^s in section III A. The operator product expansion in the Wess-Zumino- 
Witten sector can be used to replace the point-splitted product of G"'s by a sum over primary fields (z, z). Each 
(pi {z, z) in the sum corresponds to an irreducible representation present in the irreducible decomposition of the product 
representation □^•••(g)n = n- -- n®---. It is thus the primary field associated to the irreducible representation 
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with the largest quadratic Casimir eigenvalues and in the irreducible decomposition of □ (g) • • • (g) □ which 
determines the conformal weights up to subleading corrections: 



m_c^ I ^ _ (445) 
2 TV' 2 TV ^ ' 



(A^ + 2)(A^- 1)\ iV-2 



A first example is 



The conformal weights are those of the primary field transforming like the symmetric representation □□. They vanish 
for SU{2) and are negative when N > 2. Another example of a composite local operator with vanishing conformal 
weights is 0] 



(((^Lv'+V'lv'-) i^rz) [vlip+ifilifi^) (0,0)) cx (4.47a) 

^ 0. (4.47b) 



N 



N- 



k=-2N 



The confo rmal weights are those of the primary field transforming like the adjoint representation. One readily sees 
from Eq. (4.41) that the completely symmetric representation □ • • • □ — {m, • • • , 0}, which corresponds to one row 
of m boxes, has the largest quadratic Casimir eigenvalue among all irreducible representations made of m boxes (for 
example, for SU{2) it would correspond to the representation with largest possible "total spin" ^). Hence, we find 
the important result 

((^l^+)'"(z,z)(^;^_)'"(0,0)) (X (448^) 

77? N ~ \ 

h = h = - - (m^ + Nm) . (4.48b) 

As was the case along the critical line (7^ > 0, 5^ 0, there are infinitely many local composite operators with 
negative scaling dimensions along the critical line = 0, = 00. 



V. SCALING AWAY FROM CRITICALITY 



In this section, we are going to calculate how some averaged quantities such as the local density of states, generalized 
inverse participation ratios, and their spatial correlations scale with respect to energy and length. We begin with the 
dynamical scaling exponent and with the local density of states. 



A. Dynamical scaling exponent z and scaling of p{ijj) away from criticality 



The dynamical scaling exponent z relates scaling with energy to scaling with respect to length. We consider the 
local order parameter M{x) defined by 

M{x) = ('0V + w)(a;), (5.1) 

since it is coupled to the complex energy — lj ±hi in the retarded and advanced single-particle Green functions, 
respectively. Notice that M.{x) breaks the global C/(l/l) x U{1/1) invariance. We continue working in real space and 
choose a renormalization group scheme defined by the scale transformation 

X — a X, a > 0. (5-2) 

The scaling dimension of E or, for that matter, those of ijj^ (/?, A^, B^, are defined by the requirement that the effective 
action in terms of the rescaled coordinates, fields, and coupling constants has the same form as that in terms of the 
original coordinates, fields, and coupling constants. The scaling dimensions of all local fields have been calculated 
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in sections [11 and IV a long the disordered critical line > 0, = oo. The impurity strengths and are 
left unchanged by Eq. (5_^), since they couple conserved currents to vector gauge potentials. The energy scale E 
transforms like 



E = a-"" E, 



z = 1 + ^ 



1 



along the critical line 5^ > 0, gs — 00, as follows from Eqs. ( [l.43| ) and ( 3.30b ) and 



E I (fi 



h + h 



1 



9a 

TT 



(5.3) 

(5.4a) 
(5.4b) 



The dynamical exponent z allows to turn scaling with respect to length into scaling with respect to energy. When 
only Abelian vector gauge disorder is present, z = 1 + ^ is always larger than one and increases linearly with 
the impurity strength. It reaches the value z = 2 when 5^ = tt. At this value, the scaling dimension of the local 
symmetry breaking field M{x) vanishes, a signal of spontaneous symmetry breaking. The possibility of spontaneous 
symmetry breaking cannot be ruled out by the Mermin- Wagner theorem as would be the case in two-dimensional 
unitary conformal field theory, since the disordered critical theory is not unitary. On the other hand, in the presence 
of non- Abelian vector gauge disorder only, the dynamical exponent z = 2 — is always smaller than two. Hence, the 
large impurity strength limit appears to be very different depending on whether the vector gauge impurity is Abelian 
or non- Abelian. 

We now perturb the effective action Scl'^^ at the critical point > 0, = 00 with uj J dPx M{x). To lowest 
order in lo, the averaged local density of states p{ijj) scales like 



2- z 



Oico) 



(5.5) 



since, according to Eq. ( 5.4b| ), the scaling dimension of p{uj) is 2 — z. We see that in the presence of Abelian vector 
gauge-like disorder only, the averaged local density of states is finite at criticality for the impurity strength 5^ = tt. 
For larger values of 5^, p{uj) diverges in the limit lo 0. In the presence of non- Abelian vector gauge disorder only, 
the scaling exponent f3 = {2N'^ — of the averaged local density of states decreases monotonically as a function of 
N with the values P = h when N ^ 2 and /? = when N ^ 00. 



If the system is not in the thermodynamic limit but in a finite volume L , we recover Eq. (5.5) with the Ansatz 



provided 



p{uj,L) = Fp{LuL')\ujf , uj<.L- 



lim FJx) = Kp, < Kp < 00, 



(5.6) 



(5.7) 



holds. In a finite volume, the total number of energy eigenstates must be finite. This is always so (except in the limit 
00), since (3 > —1. 



B. Scaling of averaged generalized participation ratios 

In this subsection, we turn our attention to the scaling of averaged generalized inverse participation ratios for three 
reasons. First, the scaling of averaged generalized inverse participation ratios was calculated perturbatively in powers 
of e, where space is 2 4- e-dimensional, by Wegner for the non-linear a model of the metal-insulator transition at 
the mobility edge . We will obtain non-perturbatively the dominant scaling exponents of averaged generalized 

inverse participation ratios for Dirac fermions with vector potential disorder. Second, we will see that the negative 



dimensional operators found in sections III and IV control the scaling of averaged generalized inverse participation 
ratios. Thus, negative dimensional operators are not mere curiosities of the disordered critical line > Oi^s = 00. 
Finally, we reexamine arguments suggesting a relationship between the scaling exponents of averaged generalized 
inverse participation ratios and the existence of a typical multifractal wave function with vanishing eigenenergy 
(critical wave function) fl^ . 
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The concept of the inverse par ticipation ratio for the states of a disordered system was introduced to distinguish 



locaUzed from extended states |54 . Let \"^) be any state defined in the Hilbert space representing a single-particle tight- 
binding Hamiltonian. The Hamiltonian describes, say, hopping on a finite d-dimensional lattice A in the presence of a 
random site potential (Anderson model). There are |A| sites labelled byZ = l,---,|A|<oo. All states are normalized 
to one: 

leA le supp* 

Here, supp^* denotes all sites having a non- vanishing overlap with j^*). The inverse participation ratio is the functional 

Pf [*] ^ K^l*)l'' (5-9) 

le supp* 

taking values in ]0, 1]. Let now be any subset of A and define the normalized state l^'j^) by 

mn) = { ^ (5.10) 
0, otherwise. 

On the one hand, the inverse participation ratio for the state j^*^) scales with the linear size L = \ A\i as 

Pf [vl/^] = L-^^'\ r(2) = d. (5.11) 

On the other hand, the inverse participation ratio for the state where 17 has a finite number of sites in the 

thermodynamic limit L oo, does not scale with the size of the system since 

Pf - |f^r'>0, (5.12) 

is independent of L. 

(2) 

In practice, finite size scaling analysis of P\ does not distinguish extended from localized states close to the 
mobility edge. Indeed, in numerical simulations, one calculates energy eigenstates whose wave functions have at most 
isolated zeroes in contrast to, say, Moreover, for energies close to the mobility edge, the localization length 

is so large that the eigenstates, for all intent and purposes, extend throughout the systems accessible to numerical 
simulations. Close to the mobility edge, it is more appropriate to ask what is the nature of the extended state in 
the finite size system: is it a state with strongly fluctuating or uniform amplitude? Generalized inverse participation 
ratios can address this issue. Generalized inverse participation ratios are defined to be functionals on the space of 
normalized wave functions {vE'} in a volume of linear size L 



Pi'^i*] ^ / ^|vi.(,)|2., pW[vi.] ^ / (5.13) 

J supp* ^ J supp* ^ 



taking values in ]0, 1] for positive real values of g, and in [l,oo[ for negative real values of q. The usefulness of 



generalized inverse participation ratios comes from the property that P^' \^] is dominated by the largest probability 
density supj.{|^'(a;)p} when q is very positive. Conversely, when q is very negative, the smallest probability density 
inf3.{|^'(a;)p} dominates in P^ ''^'''[\E']. Thus, the q-th generalized inverse participation ratio for any extended state 
with uniform amplitude scales with the linear size of the system as 

= T{q) = d{q-l). (5.14) 

This is an example of gap scaling. ^ Any scaling departing from Eq. ( |5.14| ) signals that the state is a simple fractal 
or multifractal. This classification is based on writing 

r(g) = D{q){q-1). (5.15) 



See footnote EJ in the introduction. 
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When D{q) is independent of q, the state is said to be simple fractal, in which case gap scaling still holds. Otherwise, 
it is said to be multifractal and gap scaling breaks down. For a multifractal wave function, the supports of its maxima 
and minima scale with the system size with two different exponents D{—oo) > D{+oo). 

The functional form of T{q) is severely constrained by the fact that the wave function 5" (a;) induces a normalized 
probability measure |^(a;)p. By construction T{q) is strictly increasing and vanishes when q — 1 due to the normaliza- 
tion condition. For our purposes, the most important property of the function T{q) is that it is concave. Multifractality 
of the wave function is equivalent to the strict concavity of the function T{q). Other analytical properties of T{q) can 
be found in Ref. 0. 

In numerical experiments on some two-dimensional non-interacting tight-binding Hamiltonians with static disorder 
[ p^ , the collection of q-th generalized inverse participation ratios {Pj^^'^e v]}v calculated for some collection of 
normalized energy eigenfunctions labelled by the impurity potential V, i.e., {HQ + V)\"i' y 
that, 



P'l^E v)- It is observed 



-r(q) 



(5.16) 



for energies close to the critical energy iJ^.^, i.e., when 1 <C i <C where is the localization length. We therefore 
make the hypothesis that at criticality and in the infinite volume limit hiP'^'^^]^ ^ y\ is a selfaveraging random 
variable: 



lim 



InL 



lim 



lnp(^^[vl^^^^,^] 
-InL 



(5.17) 



(The overbar denotes averaging over the random potential V .) If Eq. (5.17) holds, it is not necessary to perform an 
averaging over the disorder to calculate T{q)^ as long as the system size L is chosen large enough. The challenge for 
a field theory describing the universal properties of the multifractal critical wave function ^E"^ v^^) ^'^ obtain the 
scaling exponents {T{q)\q £ H}. 

To illustrate the difficulties involved in calculating from a critical field theory the function T(q) defined by Eq. 
( 5.17 ), we use the definition of Wegner ||l2|] for disorder averaged generalized inverse participation ratios 



(5.18) 



and apply it to Dirac fermions with random vector potentials. Here, i labels a complete basis of energy eigenstates 
v)}i ^'^^ ^ given realization of the static random vector potential V defined in Eq (4.2). The denominator is the 
density of states per energy. It insures that V^^\lo) is normalized to one. The right-hand-side of Eq. (5.18) equals 



Pl^ [*^. v] if and only if the wave function ^E*^ vi^) been normalized to one. Although for q > Pj^ y\ 

takes values in ]0, 1[, V''''\lo) takes values in ]0,oo[, since the condition that the wave functions entering Eq. ( |5.18| ) 
are normalized has been relaxed. We now show that for energies asymptotically close to the critical energy for Dirac 
fermions, namely uo — each 'P^'^\lo) scales like |cl'|™^*-'. We will then discuss the relationship between the set of 
scaling exponents {T*{q) = zm{q)\q £ UNT} calc ulate d with the critical field theory of the previous sections and the 
scaling exponents {T{q)\q £ IN} defined by Eq. ( ^.17 ) (i.e., derived from the normalized wave functions with vanishing 
energy of Dirac fermions with random vector potentials). 

We first assume t hat, after averaging over disorder, each term in the numerator is independent of the spatial 
coordinate x in Eq. (5.18). If so, the summation over the spatial coordinates can be dropped i n the numerator, if the 
denominator is replaced by the averaged local density of states p{oj). The usefulness of Eq. ( 5.18 ) then comes from 
the identity pf 



1 



V^''Hoj)p{uj) = — lim r/«-iA(«)(a;,7]) 



q£K, 



(5.19) 



since we can relate the right-hand-side to a local operator in the critical field theory for Dirac fermions. To see this, 
recall that according to Eq. (2.4), A'^^^w.r]) is the disorder average over the g-th power of the smeared density of 
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states per volume and per energy for the complex energy E = uj + irj. ^ We now calculate the scah ng e xp onen t 
of V'^ '^^ u}) wit h resp ect to uj as we did for the averaged local density of states. We infer from Eqs. (2^), ( 3.32 ), 
( |3.34| ), ( [4.48b ), and (5.3) the transformation law obeyed by ry^^^yl(*)(0, 0) under real space rescaling on the critical 
line > 0, = oo, namely 



^{q^+Nq) q-1 ( ] Y I \ 



a > 0. 



With the help of Eqs. (5.4b), and (|5.5|), we arrive to the scaling obeyed by 7'(')(cj) to lowest order in 



w{q) = 



'(<2) 



r*{q) = D*{q){q-l), D* [q) 



9a 



N - 1 



It is remarkable that the same continuous spectrum of exponents 



w{q) 



T*{q) = D*{q){q~l), D* {q) 



(5.20) 

(5.21a) 
(5.21b) 

(5.22) 



along the critical line > 0, = 0, has been obtained by Ludwig et al. in a very different manner [pO| . They 
first found an eigenstate with vanishing energy for any realization of the Abelian vector gauge-like impurity potential. 
They then normalized the wave function in a square box of linear dimension L, thus obtaining the real and strictly 
positive wave function '^^y{x). Finally, they calculated how the disorder average |^'oy(^)P' scales with the box 
size L using the replica trick to treat the normalization constant, and assuming that this constant is self averaging. 
However, as we will discuss below, the normalization constant is not selfaveraging. Thus, the calculation in Ref. pC| ] 
cannot be viewed as a calculation of T{q) defined by Eq. (5.17). 

Using the definition for averaged generalized inverse participation ratios Eq. (5.18), we showed in Eq. ( 5.21bD that 
they obey a scaling law controlled by local operators in the critical field theory for Dirac fermions with random vector 
potentials. The scaling exponents {T*(ci)\q G IN} have been obtained non-perturbatively. By that we mean that 
T* (g) is obtained from the exact scaling dimensions of local opera tors (y ^ly,)*^ in the critical field theory (see Eq. 
( ^.35 )). Both th e sca ling exponents {T*{q)\q G IN} defined by Eq. ( 5.21b| ) and the scaling exponents {T{q)\q g N} 
defined by Eq. ( 5.17 ) characterize the statistical properties of the collection of critical wave functions {^'q yi^Y^v- 
Hence, they should be closely related although they need not be equal. In fact, they cannot be equal since Eq. ( |5.21b| ) 
implies that the scaling exponent T*{q) becomes negative for large enough q. In contrast, T( q) m ust be positive. 
There are two reasons for this obvio us difference. First, the critical wave functions entering Eq. ( ^.17 ) are normalized 
whereas they are not in Eq. ( 5.21b|) . Second, the normalization factor for unnormalized critical wave function is not 
a selfaveraging quantity for Dirac fermions with random vector potentials: 



lim I d'^x\^Qy{x)\'^ ^ lim i d^x\^Qy{x)\ 



(5.23) 



For large g, rare events, which are defined by |5'q y(a;)| being much larger than y(a;)p, dominate the average in 
p(«)(cj). Consequently, T'^'H^) 

can diverge in the limit w — > 0. 
It turns out that the is sue o f selfaveraging is far more important than the issue of normalization. Indeed, even if 
we had required in Eq. ( 5.18 ) that the wave functions be normalized to one, the corresponding scaling exponents 
{T**{q)\q G IN} could only equal the scahng exponents {r(g)|q G IN} if generahzed inverse participation ratios were 
selfaveraging. However, the parabolic shape of r* {q) strongly suggests that this is not the case, since the parabolic 
shape is characteristic of a log-normal distributed random variable ||20| , [l0| . Conversely, if the critical wave function 
is multifractal (i.e., T{q) is strictly concave), then all averaged generalized inverse participation ratios are necessarily 
non-selfaveraging. By enforcing the constraint of normalization of the wave functions in Eq. ( 5.18| ), the new set of 
scaling exponents t** [q) thus obtained must be positive, but they need not equal T{q). |^ We must then conclude 
that for Dirac fermions with random vector potentials 



^ The derivation of the combinatoric factor Cq can be found in Ref. Multiplication by r;''"^ and taking the limit r; — > 0"*" 
projects out of yl''' (tj, t;) contributions not present in 7'''''(aj). 

® The Legendre transform /**(«) = aq — T**{q), where a = has been shown to take negative values in a model of 

random resistors ISq]. On the other hand, f{a) = aq ~ T{q), where a = i^lM^ is always positive jW^. 
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{T*{q)\q e IN} ^ {r**{q)\q G N} 7^ {^('Z)I'? S IN}, 



(5.24) 



although we beheve that the strict concavity of any one of r* (q), t** (q), rjq) imply that of the two others. 

We close this subsection by reexamining an attempt made in Ref. to relate multifractal scaling exponents to 
the scaling dimensions of powers of local operator in critical field theories. Let 5*0, y(x) be a critical eigenfunction of 
the Dirac Hamiltonian with random vector potential in a square box of linear dimension L. Define the so-called 
local random event 



(5.25) 



We want to construct a normalized probability measure out of local random events. We therefore define the so-called 
local probability fO] 



Jl2 U 



Ooyix) 



We can always write 



L2 



fx 
U 



(5.26) 



(5.27) 



where e is some microscopic length scale. We made in Eq. ( ^.17 ) the assumption that \nFQy{q, j-) depends only 
weakly on L, i.e., 



lim 



InK 



InL 



InL 



g G R. 



(5.28) 



Note that Eq. ( |3.28 ) does not rule out the possibility that for any g e R averaging over randomness (denoted by an 
overbar) yields 



lim Ff^^yiq, -) = F„{q) lim - 



L 



L 



5T(q) 



6T{q) e R. 



(5.29) 



Let us make instead of Eq. (5.28) the much stronger assumption that the amplitude Fq y{q, j-) does not depend on 
the randomness V and on the system size in the thermodynamic limit: 



lim Ff 



Foiq), gGR, 



(5.30) 



i.e., /ig y{x) is selfaveraging. This is the assumption made in Eqs. (1-3) of Ref. jll]. Unfortunately, this assumption 
does not apply to our model of Dirac fermions with random vector potentials for which Eq. (pT29) with 6T{q) — 



T* (q) — T{q) holds. It is thus incorrect to use Eq. (3) of Ref. na mely to equate (see Eq. ( |5.35D ) Aq — 2 — qiA^ — 2) 
with T{q), as was done in Ref. |Q. On the other hand, Eq. ( ^.3()|) should yield a good approximation to T{q) in the 
vicinity of q = 0, since fig y{x) = exjp[q\n fiQ y{x)] ~ 1 + gln/ip y{x) is very close to being selfaveraging for small q. 



C. Scaling of averaged spatial correlations of generalized inverse participation ratios 

In the previous subsection, we have constructed a spectrum of exponents {T*{q)\q G IN} which characterizes the 
scaling of averaged generalized invers e par ticipation ratios. There exists a linear relation between the r*(<7)'s and the 
scaling dimensions {A^jg e IN}, Eq. ( 5.35| ), of some local operators in the critical field theory. In this subsection, we 



study the interplay between the operator product expansion and the strict concavity of A^ by constructing averaged 
spatial correlations of generalized inverse participation ratios. 

Following Ref. , averaged spatial correlations of generalized inverse participation ratios are defined by 



Qfa-^^)(x,y,c.) ^ ^'1^ I ' , \ ^, 91,92 >0. (5.31) 



The q-th generalized inverse participation ratio is recovered by choosing qi = q, <Z2 ~ 0- Similarly to Eq. (5.32), we 
use the identity [O 
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C. 



lim 



(5.32) 



Here, A'''^^''^^'^{x^y,uj,ri)^ which is defined by 



I Im 1 *i fTm ~l '2 

A(«-«^)(a;,y,c.,r7) = ^ —tr [G(a;, - iT?; 1/)] I j — tr [G(y, y, ^ - iry; V')] 



(5.33) 



should be compared to Eq. {2A). We calculate the scaling of yl^'^i '*2)(a;, y, w, 77) with respect to /arye separation |x — ?/| 
and with respect to small to ^ Q. In other words, we work in the regime 



e < |a; - y| < L^, cx \uj\~ 



(5.34) 



Here, e is the ultraviolet cutoff, say, the lattice spacing. Of course, 77 > can be taken as small as needed. Scaling 
is calculated as follows The first step is to perform a rescaling which brings the separation down to the lattice 
spacing e. The second step makes use of the operator product expansion. The final step is to perform a rescaling 
to energy scales uj and fj which are arbitrarily chosen. Under each rescaling, we use the scaling dimensions of the 
appropriate operators. In particular, to lowest order in w, we use the scaling dimensions at criticality (see Eqs. 
( |3.30b| , |4.48b| )) 



(5.35) 



for the operators {if\_(p_\ [x), lip\.ip_j (y), and lip'^ip_\ (e), respectively. 
Under rescaling, A^'^i^''^) (a;, a;, 77) transforms like 

,9l , , / + \<l2 



according to section HI. Hence, under the rescaling 



(5.36) 



ax, y 



ay, uj 



a uj, rj ^ a rj, a 



\x - y\ 



«1, 



(5.37) 



we find that 



A('?i-«2)(a;,y,c^,77) 



k - y\ 



\x - y\ 



A +A 



9l+«2 



(y) 



(5.38) 



On the second line of Eq. ( ^.3^ ), we have performed an operator product expansion from which only the ttios^ 
relevant operator of the expansion is kept. Notice that the identity operator is not the most relevant operator in the 
operator product expansion. It would have been, had we considered, say, the pure Gaussian field theory with the 
hierarchy of operators exp( gi0) a nd the positive and convex spectrum of exponents {q^\q G IN}. The two factors on 
the right-hand-side of Eq. ( 5.3^ ) both depends on the separation |a; — y\. However, the rescaling to the arbitrarily 
chosen energy Q < ui 



by, Lo ^ b ^uj, rj ~ b ^7), b 



(5.39) 



yields 



A^''^^'^^\x,y,u,Ti) 



A +A -A , 

e \ '1 «2 ''i+''2 



\x - y\ 



A +A -A , 
£ \ '1 '2 ''l+''2 



F ~ y\ 



\ \ / / r,(1.0),- .-s 



i-A„ 



(5.40) 
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The scahng function F^'^'- {u) , fj) is defined by Eq. (5.40). It is giv en by the expectation value at energy (D — i?7 of 
the most relevant operator in the operator product expansion in Eq. ( 5.38 ). By assumption, it does not depend on y. 
For given qi2, it therefore only depends on uj and fj. Calculating F^''^'^^^{uj,fj) is beyond the scope of our approach. 
However, some of its limiting behaviours can be inferred from the fact that ''^2) (cj, fj) enters an averaged generalized 
inverse participation ratio. 

We consider the limit a; > 0, ^ 0. Because we are working in the thermodynamic limit at fi nite real energy w, 
the energy eigenstates must be localized. Hence, Q(''i'*2)(2;, oj) must be finite so that from Eq. ( f3.32 ) 



91+92-1 



This gives us the scaling 

Q('^-«2)(x,y,c.) Pico) 



F - y\ 



A„ +A„ -A„ 



^,1+<J2 +^(91+92-1) 



(5.41) 



(5.42) 



This scaling has two important features. First, the exponent characterizing the scaling with separation is always 
positive since 



A -A ^ 

92 9i+92 



Secondly, the exponent characterizing the scaling with respect to lo is that of the {qi 
participation ratio: 



(5.43) 



(72)-th generalized inverse 



1 r 

z 



h<32 



z {Qi +92-1) = ^(9i + 92) + 



2- z 



(5.44) 



Consequently, we recover the scaling for generalized inverse participation ratios when either or vanishes as it 
should be. 

Eq. ( ^.42 ) tells us that the property of strict concavity A^^ + A^^ — A^^^^^ > 0, (7^,(72 > 0, insures that spatial 
correlations decrease with separation. Negative dimensional scaling exponents also control the scaling with respect to 
energy, or, equivalently, with respect to the effective system size. This latter scaling implies that the normalization of 
the local observable can blow up with the system size in the thermodynamic limit. In fact, Eq. ( ^.42 ) reproduces the 
scaling for spatial correlations of so-called non-normalizable box-observables |3^,|ll|,|l^ . Also note that the identity 
operator never enters as the most relevant operator in the operator product expansion between pairs of local operators 
with scaling dimensions A^^ and A^^ , and this, in turn, is guaranteed by A^ becoming negative for sufficiently large 
q. In contrast, it is much more difiicult to find a critical point characterized by an infinite set of positive scaling 

Mn > is generated from x„ and x„ as 

^1 '^2 9i ^2 



exponents {x > 0\q G IN}, which are not Hnearly related, and such that 



the most relevant operator in the operator product expansion of their associated operators [ pl] |. 

In summary, several scaling properties of the multifractal measure induced by the critical wave function of Dirac 
fermions with random vector potentials are shared by the probability distribution corresponding to averaged general- 
ized inverse participation ratios. In the literature it is assumed that the strict concavity of T*{q) for small q implies 
that of T{q). Although this assumption seems reasonable, it would be highly desirable to verify it by reconstructing 
{T{q)\q G IN} from {t* {q)\q G IN}. Moreover, by doing so it could be established whether or not the existence of 
negative dimensional local operators implies multifractality of the critical wave function. 



VI. CURRENT ALGEBRA FOR GRADED SYMMETRIES: (7(1/1) CURRENT ALGEBRA 

In the previous sections, we have studied some disordered critical points through their fixed point Lagrangians. In 
this section, we take a very different approach. We want to describe new disordered critical points using conformal 
and internal symmetries and the associated current algebra. We will see how the current algebra approach allows 
us to access new fixed points which cannot be reached perturbatively from some known fixed point (contrast this 
approach to that in Refs. ]l8| , p^ ). 

Current algebra (and the associated bosonization technique) is a powerful method for studying strongly correlated 
electrons in one spatial dimension | ]5^ . Current algebra is also very useful for studying critical points in two dimensions 
which are characterized by certain continuous symmetries. Effective theories for non-interacting disordered systems 
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have a graded symmetry (or supersymmetry) [ p8| . For example, we have seen in section || that disorder average over 
products of Green functions can equally well be represented by path integrals over fermionic or bosonic coherent states. 
In this section, we study the current algebra of graded symmetry in two dimensions. We classify all critical points 
for two-dimensional systems with unbroken [/(l/l) graded symmetry. We recover the critical points corresponding 
to U{1) vector gauge-like disorder studied in the previous sections. We also find new critical points. The operator 
content for these critical points is constructed. In particular, we show that most critical points arc characterized by 
infinitely many local operators with negative dimensions. 

In the current algebra approach, one first constructs the commutator and anticommutator algebra for the conserved 
currents. At a critical point, conformal symmetry allows one to calculate all correlations between the currents from 
the commutator and anticommutator algebra. In this way, the theory is completely specified. 

A. Current algebra description of a free model 

We begin with the simplest model with graded [7(1/ 1) symmetry. It is given by the partition function 



£0 = ^ + ifi+dgip^ + ip^_dzip_ + (p^_^dzip_^ . (6.1b) 

Here, (x,r) denotes the coordinates of two-dimensional Euclidean space, and the fields are the chiral components of 



Dirac spinors introduced in Eq. (3.15|). For convenience, we have rescaled the spinors by the factor vStt. The only 



non-vanishing two-point Green functions are 

(Vl(^)V'+(0)> = {^+{z)^lm - i, (6.2a) 

-{^l{z)ip+m = - ^, (6.2b) 

(V'i(^)^_(0)) = (V'_(^)VA(0)) = i, (6.2c) 

z 

-{^Uz)^_{0)) ^ {^_{z)^U0)) = i (6.2d) 

z 

where z = x + iT and z = x — it are the independent holomorphic and antiholomorphic coordinates, respectively. 

1. Symmetries of Cq 



The action in Eq. (6J_) is left invariant by the global chiral transformation 

- e--< ^'l V+ = e+'< (6.3) 
where 9^ is any arbitrary real number. The associated current 

Jt{z) = (V'U+)(z), (6.4) 
is conserved, i.e., its conservation is equivalent to it being a holomorphic function 

d, Jtiz) = 0. (6.5) 
Similarly, there are three additional independent global U{\) chiral symmetries with the conserved currents 

jt{-z) = (vi^-)(^), jn^) = jn^) = {^iv-m. (6.6) 



respectively. The action in Eq. (6.1) is also left invariant by the four independent infinitesimal transformations 



5ri'P± = [S±l±,i^±] = +V±d±, \V>± = [S±1±.V±\ = -0±^±, (6.7a) 
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and 

5^4 = [0±V^,vi] = -4^±> = [0±V±,i^±] = +0±V±, (6.7b) 

respectively. The four infinitesimal Grassmann numbers 9_, 0_^, and are independent. The associated 
conserved fcrmionic currents are 



= V'+V'+, V- = 'y'iV'-, V+ = '4^+'P+, V+ = iplf-- (6.7c) 



Taken together, these eight conserved currents, of which four are bosonic and four are fermionic, generate a graded 
?7(1/1) X U{1/1) symmetry. The bosonic generators of the diagonal subgroup U{1/1) are the currents and with 



2{j, + uJ - \J-J - \jt + Jt)' 2\], + ijJ - \j-J - \jt-Jt 

The fermionic generators of the diagonal subgroup [/(l/l) are the currents rj^ and 77^ with 











(::) 













(6.8) 



(6.9) 



Each component of the bosonic and fermionic currents generating i7(l/l) is separately conserved due to the chiral 
invariance. 



2. Symmetries of Zg 

The global C/(l/l) x C/(l/l) symmetry of Cq is shared by the measure in the partition function. Hence, Zq is also 
left invariant by global [/(l/l) x [/(l/l) transformations. The correlation function {A^_^_ ■ ■ ■ A.^-^^B^_ ■ ■ ■ -B„_), where, 

say, are defined in terms of V'+iV^+i </5+,</5+, factorizes according to 

{A,^ ■ ■ ■ A^+B,_ ■ ■ ■ B„_) = {A,+ ---A^+){B,_---B^_), (6.10) 

due to the decoupling of the + and the — sectors. Moreover, the correlation function corresponding to the + subscript 
must be holomorphic, whereas the correlation function corresponding to the — subscript must be antiholomorphic. 
Without loss of generality, we will only consider the holomorphic sector of the theory. The subscript + will be 
suppressed from now on. 



3. Operator product expansion 



We now study the operator product algebra with the help of the operator product expansion (OPE) 
two operators, A{z) and B{z). The OPE of the product A{zi) B{z2) is the Laurent series expansion 



Consider 



A(zi) B{z2) 



{zi - Z2)' 



:0n{z2) + ••• + 



Zl - 22 



-01(22) + 0{1) 



(6.11) 



Only the poles in the limit of zi — > 22 are kept in the OPE. In the following, we will always omit the non-diverging 
term 0(1) in the OPE. For example, the operators tj}\ ■ ■ ■ ,ip have the following OPE: 



iIj\z)iP{0) = ^(2)V'^(0) = Az)^HO) = -'P\z)ip{0) 



(6.12) 



Another example of an OPE is that between the currents (J, j, r], ff) generating the [/(l/l) symmetry. However before 
calculating them, the quantum nature of the currents must be dealt with care. In particular, we must remove any 
divergences which arise when multiplying operators at the same point. In the context of the cur rents, such divergences 
are removed by subtracting from the naive definition of the classical currents in Eqs, 
between the ■ ■ ■ ,(p. For example, the precise definition of j{z) is not given by Eq. (|6.8 



|,|J) the poles in the OPE 
, but by 
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lim 



'4)\z + e)'4}{z) — ifi^ {z + e)ip{z) 



■.^P\z)^{z) 



ip\z)Lp(z) : 



(6.13) 



The normal ordering : • • • : extracts the finite contribution from the operator product defining the current j , since 
one can show that the Hmit e — > in the correlation function 



[il;\z + e)i;{z) - ^\z + e)^{z)\\{0,{zA - -([]o,(z,) 



(6.14) 



is finite and thus defines the correlation function {j{z)Y\iOi{zi)) . With this definition of the currents and with the 
help of the OPE between • • • , (/j}, we find the OPE between the currents to be given by 

(6.15a) 
(6.15b) 
(6.15c) 
(6.15d) 

An alternative definition of the OPE for the current s can be found in appendix [A] 

The operator algebra specified by the OPE of Eq. ( ^.15| ) is very powerful. Indeed, from the OPE alone, i.e., without 
any reference to Zq, we can calculate all the correlation functions between the currents (J, j, 77, 77). Here, we assume 
that expectation values of all currents vanish 



J(z)J(0) = 


0, 


^(^)j(0) = 


\, J(z)?7(0) = 0, J(z)77(0)=0, 




i(^)^(o) = 


2 

^' 


mm - 


0, ]{z)7i{Q) = -\(0), ]{z)n{Q) = 


2 , . 
+ -^(0), 


r;(z)J(0) = 


0, 


mm = 


+ -?7(0), r7(^)'7(0) = 0, 77(2)77(0) = - 

z 


Z'^ z 


n{z)J{Q) = 


0, 


n{z)]{o) = 


--77(0), ^(z),7(0) = +1 + i J(0), 
z z^ z 


mm) - 0. 



(J(z)) = {j{z)) = {^{z)) = {fj{z)) = 0. 



(6.16) 



In this sense, the OPE of Eq. ( |6.15 ) completely specifies and defines our model. To illustrate the power of the OPE, 
let us consider the correlation function {j{zi)fj(z2)'ri{z3)). We know that this three-point function is a holomorphic 
function of z^'s. We may view it as a function of zi. If so, this function must have poles at Z2 and Z3 and must vanish 
asymptotically as zi — > 00. The poles are then completely fixed by the OPE: the three-point function has a first 
order pole at Z2 with residue 2(77(^2)77(^3)), and a first order pole at Z3 with residue —2(^(22)77(23)). This completely 
determines 



(^(21)77(22)77(23)) = 
With the help of (^(22)7/(23)) 



-(77(22)77(23)) 



-(77(22)77(23)). 



21 - 22 - Zl- Z3 

2-, we find that the three-point function is given by 

2 



(22-23) 
(^■(21)77(22)7/(23)) 



(6.17) 



(6.18) 



(21 - 22) (21 - 23) (22 - 23) 

It is important to realize that not every possible OPE gives rise to a consistent theory. For example, to calculate the 



three-point function in Eq. (6.18), we expressed it in terms of two-point functions whose coefficients are determined 
by the OPE between the currents. But by so doing, we chose a specific decomposition path corresponding to viewing 
the three-point function as a function of 21. One needs to show that different decomposition paths always lead to 
the same correlation function. In the above example we may instead first view (j(2i)77(22)77(23)) as a function of 22, 
which has a first order pole at 21 with residue —2(^(21)7/(23)), a second order pole at 23 with residue (j(2i)), and a 
first order pole at 23 with residue (j(2i) 7(23)). Thus, 

2 



{j{zi)7j{z2)lliz3)) 



Z2 - Zl 



{m)m)) + 



Z2 - Z-i 



{3{zi)J{z3))- 



(6.19) 



Comparing Eq. ( |6.17D with Eq. ( |6.19| ), we see that the three-point function can be decomposed into different 
Laurent series (OPE). Thus, the OPE between the currents need to satisfy certain self-consistent (associativity) 
conditions insuring that all current correlation functions are uniquely defined. In our example, this is the case. The 
self-consistency of the OPE in Eq. ( |6.15 ) follows from it being derived from the partition function Zq. 
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B. General (7(1/1) current algebra for interacting models 



We are now ready to describe the philosophy behind the current algebra approach. We know that the model specified 
by Zq in Eq. (^^) describes a critical point characterized by algebraic decaying correlation functions. Moreover, the 
representation of Zq in terms of the free fields ■0+ j ■ ' • ; 'P- allows us to easily calculate all properties of the critical 
point. On the other hand, for a generic interacting theory such as 



r — 



/ D[^l^^] e^-^'^'^ (6.20a) 
ft + (VA0_ + ^l^^) + (vU- - + • • • , (6.20b) 



the interactions drive the system to a new fixed point. The interaction with coupling constant > can be thought 
of as being induced when integrating over Abelian vector gauge-like disorder in the effective action S'^^'^^ of Eq. (3.25). 



Similarly, the interaction with coupling constant g^^^ > corresponds to integrating over a randomly distributed mass 
for Dirac fcrmions. However, from the point of view of the current algebra, gj^ can be negative as we will see 
later on. Finally, "• • •" represents other interaction terms. We would like to know what are the critical properties 
of the new fixed point. The difficulty is that the fields tA+i ■ ■ ■ j V- do not provide a convenient description of the 
new fixed point, making it extremely difficult to calculate its critical properties in terms of • • • , <^_- However, it 
is sometimes possible to take advantage of the fact that the interaction C-^^^ has a f/(l/l) graded symmetry. At a 
critical point of the interacting model, the C/(l/l) symmetry is promoted to the U{\/1) x C/(l/l) symmetry in the 
low energy sector, due to the decoupling of the holomorphic and antiholomorphic sectors at low energy. In this case, 
many critical properties of the new fixed point can then be derived solely from the knowledge of U{\/\) x C/(l/l) 
current algebra. The current algebra approach to a critical point in two dimensions is based on two assumptions: 
i) the critical point has a conformal symmetry (which is responsible for promoting the [/(l/l) to U{1/1) x [/(l/l) 
symmetry), ii) the conserved currents have algebraic decaying correlations, i.e., they belong to gapless sectors. 

To construct the critical points with U{\/\) x U{1/1) symmetry, we denote the infinitesimal generators of the 
U{l/l) graded symmetry in the holomorphic sector by {1-^,12,1^, I /^){z) = {J, j,r],f]){z). Without loss of generality, 
we only work in the holomorphic sector. Since conserved currents cannot develop anomalous dimensions, the conformal 
weights of /j,, a = 1, • • • , 4 are (1,0). Consequently, the leading pole in the OPE among the {I^} must be of second 
order, and the most general OPE between the holomorphic components of the generators of the f/(l/l) graded algebra 
takes the form 

4(z)4(0) = ^ + + 0{1), a,6=l,...,4. (6.21) 

The coefficients k^^^ of the second order pole contain information about the dynamics of the critical point. We will 
see below that they are severely restricted by the condition that fc^j, be an invariant second rank tensor with respect 
to the U{1/1) symmetry. The coefficients /^j,^ are the structure constants of the graded Lie algebra U{1/1). 
To see this last point, we introduce the conserved charges 



1 



Qa{T) = — dxIaix,T), a = 1, • • • , 4, (6.22) 

which generate finite J7(l/1) transformations. Recalling that correlation functions defined by a path integral corre- 
sponds to r-ordered correlation functions in the operator formalism psf , we find (a, 6 = 1, • • • , 4) 

^[Qa(0),Qb(0)]±nO(a;.,T.)^ = (^^) l^j dzjjz'I^{z')I,{z)\{0{zi)^, (6.23) 

where the integration contou r C i s a small circle centered at z. We have used the fact that is a holomorphic function 
of z in the last line of Eq. ( |6.23| ). Care has to be taken in Eq. ( |6.23 ) to use the anticommutator [•,•]_!_ whenever 



and are both fermioni c, wh ereas for al l oth er cases the commutator [ • , • ]_ is to be used. Substitution of the 



OPE between currents, Eq. (3.21), into Eq. (6.23) yields the desired result 



[Qa.Qb]±^ labcQc. a,6=l,---,4. (6.24) 

We see that the first order poles in the OPE of the conserved current are determined solely by the symmetry group, 
and thu s do not depend on the interaction. For the J7(l/1) symmetry, can be read off from the first order poles 



in Eq. (3.15). Moreover, we can generalize the calculation leading to Eq. (6.24) by introducing 
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1 



+ 00 



Q^(t) - — y dx Ux,t) f{z), a=l,...,4, (6.25) 

for an arbitrary holomorhic function f{z). We find 

[Qi{T),I,{x,T)]^ = r{z)k,, + f{z)UJ,{z), a,6=l,-.-,4, (6.26a) 

or equivalently 

[I,{x',r)J,{x,T)]^ = -S'{x'-x)k^, + S{x'-x)UMz), a.b^l,---^- (6.26b) 

To be consistent with the [/(l/l) symmetry, fc^j, must be an invariant tensor of the symmetry group. This constraint 
alone requires fc^^ to depend on only two (real) parameters (fc,fc^), which are uniquely determined by the (7(1/1) 
symmetric interaction. To see this, we consider first an infinitesimal transformation generated by Or] + 6fj. Under this 
transformation, the currents (J, j, 77, fj) change by the infinitesimal amounts 

{S,, + S^)J = 0, {S^ + S^)j = 2fje + 29r], (6^+6^)7^ = 0J, {S,^ + S,^)f] = 9. (6.27) 

Here, as a result of current conservation, the transformation law ob eyed by the currents is independent of the inter- 
action, and thus can be derived from the free theory Zq of Eq. (6.1). Now, by requiring that fc^^ be invariant under 
U{1/1), it follows that any infinitesimal chan ge S{I „{z)I^{0)) generated by the currents does not contain poles higher 
than first order in z. For example, from Eq. (6.27) one finds 



{S^ + S^){rj{z)jiO)) - ej{z)j{0) + 2rj{z){m0 + HO)) 

0{-), (6.28) 



' {k jj + 2kjjfj) 26'fc,pj 



z 

which implies that kj^ + 2k^f^ ~ and k^^ = 0. Repeating this argument for all S(^I^{z)If^{Q)) , we find that the most 
general OPE between the currents generating [/(l/l) transformations in the holomorphic sector are given by 

2k 

J(z) J(0) = 0, J(z)j(0) = J(z)r](0) ^ 0, J(z)fy(0) = 0, (6.29a) 

z'' 

2k k ■ 2 2 

j{z)JiO) = —, j{z)j{0)^^, j(2h(0) = — ^(0), jiz)fjiO)^+-m, (6.29b) 

Z'^ Z'^ z z 

r,(z)J(0) = 0, r7(^)j(0) = +-?7(0), 77(2)77(0) = 0, 77(^)^(0) = + --^(0), (6.29c) 

Z Z'^ z 

^(z)J(0) = 0, 77(z)j(0) - --77(0), 77(2)77(0) = +4 + -^(0), 77(2)77(0)^0. (6.29d) 

Z Z'^ z 

We show in appendix O that A: = 1, /c - > corresponds to the current algebra derived from the effective partition 



function zir'°^ of Eq. (|3.25D with C/(l) vector gauge-like disorder. In other words, the [/(l/l) x [/(l/l) current 
algebra given by Eq. ( |6.29p with fc = 1,A:^ > describes a two-dimensional random [/(I) vector potential model at 
criticality. It is also shown that the OPE between the currents is self-consistent (associative) for all real values of 



(fc, fcj). In other words, all correlation functions for the currents are uniquely defined by their OPE in Eq. (6.29). 

We close this subsection with the construction of the energy-momentum tensor associated with the OPE in Eq. 
( |6.29 ). The holomorphic component T of the energy- momentum tensor (EM) describes important properties of the 
fixed point. We sketch here how the energy-momentum tensor T in the holomorphic sector is uniquely specified by 
the requirements that: i) T is Hermitean and is solely constructed from the currents ii) T has conformal weights 
(2, 0), in) T transforms like a singlet under [/(l/l). At the classical level, conditions i) and Hi) are met by the Ansatz 

4 

T{z) = {laTabh){-A. (6.30) 

a, 6=1 

where T^^ is Hermitean. The condition that T be a [/(l/l) singlet is purely classical. Imposing Hi) yields 

T{z) = K^{Jj + 7777 - 7777) (z) + K^{JJ){z), (6.31) 

where 2 ^■re real. At the quantum level, conditions i) and ii) must be treated with care. The quantum definition 
of T involves taking the limit 
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T{z) = lim <^ Ki 



i J(z + e)j(z) + ij(z + e)J{z) + r]{z + e)ri{z) - r]{z + e)r]{z) 



+ K2 J(2; + e) J(z) 



(6.32) 



which defines the product of operators at the same point. In general, we need to subtract out some infinities, i.e., the 
poles that appear in the OPE o f the currents, in order for the E M ten sor to be finite in the limit e — > 0. However, one 
can show from the OPE in Eq. ( 6.2£ ) that T, as defined b y Eg. (|6.32| ), is already finite so that additional subtractions 
are not needed. Also note that the right-hand side of Eq. (6.32) is a C/(l/l) singlet even for finite values of e. The two 
additional coefficients k.^ and can be determined from the OPE of the EM tensor. In a two-dimensional conformal 
field theory Go|, the OPE of the EM tensor must have the form 



T(z)T(0) = ^ + 4t(0) + -d,T{Q). 



(6.33) 



The residue of the fourth order pole determines the so-called Virasoro central charge c. The residue of the second order 
pole determines the conformal weight of the energy momentum tensor, namely two. The constraint ii) is equivalent 
to imposing Eq. ( ^.3S ). One can show from the OPE in Eq. (6.29) and the definition Eq. (6.32) that Eq. ( 6.33| ) 
holds if 



1 

2fc' 



4- fc, 



(6.34) 



in which case the Virasoro central charge vanishes: c — Q. One can also show that the currents {J,j,'ri,'i]) have 
conformal weights (1,0), since their OPE with T is given by: 



1 



1 



T(z)4(o) = -^4(0) + -djM 



l,---,4. 



(6.35) 



Conversely, we show in appendix^ that, requiring Eq (6.35) to hold together with the definition Eq (6.32), uniquely 
defines the energy momentum tensor. 



C. Operators in one- and two-dimensional representations of the f/(l/l) Lie algebra 

The U{l/l) current algebra contains many operators apart from the four currents {J,j,ri,fi). For example, from 
the fact that the current J is a singlet of [/(l/l), we can use J to construct operators which form a one-dimensional 
representation of the [/(l/l) Lie algebra. Let us consider the following operator 

f^{z) EE e-'*'^'^^-\ (6.36a) 

where q is real, and 0j{z) is the solution of 

d,9j{z) = iJ{z). (6.36b) 

We use the OPE between the currents to construct the [/(l/l) charges of and its (holomorphic) conformal weight 
hj: . One verifies that /^(z) has the following OPE with the currents 

J(z)/,(0) = 0, j{z)f^{0) - -^/,(0), v{z)fg{0) - ^(^)/,(0) = 0. (6.37) 



The OPE in Eq. ( |07| ) gives the [/(l/l) charges {J,j,r^,ff) = (0,-g,0,0) of Clearly, /t/^ is then a [/(l/l) 
singlet, and carries a one-dimensional representation of [/(l/l). The OPE in Eq. ( |6.37 ) allows us to reduce any 
correlations between currents and /^'s into one that involve only /^'s. The correlation between /^'s is given by 

/-g(^2)> = constant, (6.38) 

since J(z)J(0)=0. This implies that has vanishing conformal weight: hj = 0. More generally, the n-point 
correlation function for is given by 
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114 (^^) 



ifEr=i9..y^0, 
constant if 1i — 0- 



(6.39) 



The U{l/\) current algebra also contains operators which form two-dimensional representations of the [/(l/l) 
algebra. To see this, we define 



(6.40) 



wherep is a positive integer. By showing that tpj^ijjp + (pjj(pp is a singlet under ?7(1/1), it follows that {ipp, (pp) transforms 
like a doublet of [/ (l/l ). Notice that this is certainly true when p = 1. The J and j charges of {ipj,, (pj^, tpp, ipp} can 
be read from Eq. ( |6.40| ). They are {p,p, —p, —p} and {2 — p, —p,p — 2,j5}, respectively. Hence, on the one hand, the 
J and j charges of V'pV'p ^ind fpfp both vanish. On the other hand, because of the infinitesimal transformation law 



(6.41a) 
(6.41b) 



which is generated by the fermionic currents t], and fj of C/(l/l) (compare with Eq. ( |6.7| )), only the sum of V'pV'p ^-nd 
ipjjfp is a C/(l/l) singlet. Moreover, we can read from Eq. ( |6.4l|) and the J and j charges of '/'p, V'pi fp^ the OPE 
between the currents and {-0^, (^J,, ■(/'p, ipp}. They are 

^0t(o), ^(z)v4(0) = ^^^(0), ^(z)Vj(O) = 0, 
J(z)^t(o) ^ P t (0), j(^)¥'^(0) = -^yp{0), viMi^) = 0' '7(^)¥'U0) = 



J(z)^t(o) ^ ^^t(o), j(^)^t(o) 



J(z)0p(O) = -^0p(O), j(^)^p(0) = -^V^p(O), r?(z)0p(O) = 0, r7(z)Vp(0) = ^¥>p(0), 
J(z)^p(0) = ^^^iPpiO), j{z)ipp{0) = ^^p(0), r7(^)¥'p(0) = -^V'p(O), j](^)¥'p(0) = 0. 



(6.42) 



We see that the OPE of (3.37) and (3.38) are well defined even when p is not an integer. From now on we will regard 
p to be a positive real number. 

Even more general doublet can be constructed if we combine the ■ ■ ■ , (pp with fq by defining 



^Iq = fl^l^ "flq = /gVp, 'i'pq = ^pfq, Vpq = ^pfq- 



(6.43) 



Here, p is a positive real number while q is real. The OPE between (ippqjPpq) and the currents are very similar to 
those between ('0p,(^p) and the currents, since the multiplication of fq merely shifts the quantum number of j. In 

fact, the OPE between {ippq,(ppq) and the currents are still given by Eq^ (S.42) except for p in the OPE with the 

current j being shifted to p — q. In fact, the transformation law of Eq. ( 3.41) also applies to ( •i/j^p, (f^^). 

From the OPE of Eqs. (6. 2£ ,6.37, 3.42), together with the definition of the EM tensor in Eq. ( |6.32D , we can calculate 
the OPE between T and {ippq, fpq)- It is 



T(z)^p^(0) = -^%q{Q) + - (9>p,) (0), 
r(z)^p^(0)=.%^p^(0) + i(9,^p,) (0). 



Eq. (6.44) implies that {ippq, ^Ppq) are primary fields [Q with (holomorphic) conformal weight given by 



P_ 
2k 



(6.44a) 
(6.44b) 

(6.45) 



For the free theory of Eq. (3J), {k,k-) = (1,0). In this case Eq. ( |6.45| ) predicts /i^q — f' which agrees with the 
dimension of ^jjipP~^ and ip^ calculated using the dimension h^p ^ = ^. Similarly, (V'pqiVpq) have the same scaling 
dimension hpq. This is a consequence of the charge conjugation symmetry of the theory. 
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We have shown that the OPE between the currents in Eq. (6.29) completely determine all their n-point correlation 
functions. We want to know if the same holds for the doublets (V'pq, Vp?); whether it is possible to calculate any 
n-point correlation function between those new operators and the currents. We begin with the two-point correlation 
functions between , (/^p, , V'pq j Vpq • On the one hand, from Eq. ( |6.42[) , V'pqV'pg and fpqflq have vanishing (J, j) 
char ges, a nd, therefore, can acquire non-vanishing expectation values. We see that current conservation together with 
Eq. ( 6.44 ) requires that 



(6.46) 



On the other hand, Eq. ( 6.42 ) tells us that i'pq'-Ppq and 'Ppqtplq have the (J, j) charges (0, —2) and (0, +2), respectively. 
Hence, they cannot acquire non-vanishing expectation values due to current conservation: 



This allows us to relate the real parameters and C^. Indeed, if we use the transformation law 

Sf^i^lq = 0, S,-^iplg = -^ijlg9, 5,-^^pq = +\/pdll)pq, Sfjifipq = 0, 

under infinitesimal transformations generated by 77, and fj, we can show that 

1 



(6.47) 

(6.48a) 
(6.48b) 



{^ + S,){^pq{z)^PlqiO)) = ^9 i-C^+C^) 



(6.49) 



Eqs. ( |6^ ) a nd (|6.49| ) then imply that = C^. Setting C. 
reduces Eq. ( |6.4(: ) to 



1 through rescaling of the doublet (?/> , if. 



i^pqi^Hqi^)) = Mz)^lm - 



pq 1 rpqj 1 



(6.50) 



It is useful to consider the OPE of two doublets {ippiqn ^piqi) and ('0p2«2i '/'p2<?2)- One can easily see that the 

product {lppiqi,iPpiqi) X (V'p^ 

92 7 fp2q2) contains a new doublet {il'paqa i Vpaqs) with — pi+p2 and 93 — 91 + (72- Thus, 
the OPE of (V'pigi J '/'piiji ) and (V'p2 92 7 '/'p2<?2) niust be of the form 



<Ppi9l(^)V'p292(0) 
V'pigi(^)¥'p292(0) 

fpiqi (^)V'P2<Z2(0) 



1 



zhi+h2-h:,'^P3l3 H : 

^hi+h2-h3'^P3l3 ^ ■ 



(6.51a) 
(6.51b) 
(6.51c) 



where hi are dimensions of (ippiqij 'fpiqi)- This formula will be used below to discuss the locality of operators. 



D. Local operators and quantization of k 

According to the general principles of two-dimensional conformal field theory (GET) ||4C| ], any generic operator can 
be written as the product of operators in the holomorphic (z) and antiholomorphic (z) sectors: 

0(z,z)=0,(z)0,(z). (6.52) 

The OPE of such mixed operators are given by 

O,(z,z)O,(0,0) = ,h^-h.-h, ^h.-h.-K 0^(0,0), (6.53) 

c 

where (/iq, ha) are the conformal weights of Oa, and C^^ are complex numbers. 

The concept of a local operator is very important in our theory. In CFT, two operators Oa(z,z) and Ob(z,z) are 
said to be mutually local if their OPE (i.e., the right-hand side of Eq. (6.53)) is single- valued. This implies that the 
spin difference 
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Sc-Sa-Sh= {K - K) - {K - K) - {K - K) (6-54) 

is an integer for all c with non- vanishing C^^. A set of operators is said to be made of local operators if they are all 
pairwise mutually local. We will call the set of all local operators the "center'^ of the operator product algebra. From 
this definition we see that the correlation functions between local operators (i.e., between operators in the center of 
the operator product algebra) are always single- valued functions. 

With these definitions in hand, we come back to the theory defined by Eq. ( 6.2C| ), assuming only that the interaction 



be invariant with respect to U{1/1). In the presence of a generic U{1/1) symmetric interaction c\^Ji^^\ the spatial 
dependency of the doublets (■0+,(^+) and (■0_,(,c_) are no longer purely holomorphic and antiholomorphic, re spec tively. 



This is not to say that the holomorphic and antiholomorphic sectors are absent. Any critical point of Eq. (6.20) with 
U{1/1) symmetry still possesses two decoupled sectors characterized by holomorphic and antiholomorphic correlation 
functions. Rather, the doublets {'ipj.,ip ±) are not necessarily primary fields anymore. Instead, they are of the mixed 



form in Eq. (6.52). (They do, however, still transform like C/(l/l) doublets due to current conservation.) 



It should be obvious that i/)^ , • • • , and their products in our model Eq. ( |6.20D arc local operators provided 
we assume that the interaction C'^J^^^\ besides being [/(l/l) symmetric, is also local (i.e., only involves product of 
fields at the same point). If so, the doublets {ip^, (pj.) and (?/'±, (p±) are local regardless of whether the model is at a 
critical point or away from a critical point and regardless of the nature of the critical point as long as it possesses the 
C/(l/l) symmetry. Thus, any critical point associated to C^J^^^"^ must contain a set of local operators that carry the 
same U{1/1) quantum numbers as tj}±^'P±^''P±')'P± ^^'^ their products do. We recall that C/(l/l) denotes the diagonal 
subgroup of U{1/1) X U{1/1), which appears at the critical point. This puts a strong constraint on the possible critical 
points that ^Cj^j^^^^-* can have. We now show that the constraints of locality and of current conservation restricts the 
possible values of the central charge k to be the inverse of some integer m: 

fc = — , me Z. (6.55) 
m 

We also show that all the critical points with [/(l/l) x [/(l/l) symmetry for which ■ ■ ■ are local are labelled 
by two integers m S Z and one real number kj. 
The proof uses the components 

Jt = HJz- Jz) : Jr = i U z " J z) > (6-56) 

which were defined in Eq. (|6^). Here, and are the bosonic generators of the diagonal subgroup [/(l/l) 
of [/(l/l) X [/(l/l). For convenience, we have listed some [/(l/l) x [/(l/l) charges in table |[ from which one 
immediately finds that the charges of (V'+, "0-, <P_) are i(— 1, —1, 1, 1), while their charges are i(— 1, 1, 1, — 1). 
Let us assume that we are at a critical point with [/(l/l) x [/(l/l) symmetry characterized by the real valued Kac- 
Moody central charges k and kj . We know from the previous subsection that there exists primary fields /^^ and /^^ 
{q real) with vanishing conformal weights and carrying a one-dimensional representation of the [/(l/l) Kac-Moody 
algebra in the holomorphic (z) and antiholomorphic (z) sectors, respectively. There also exists primary fields tppqz, ^pqz 
and i^pqzj'-Ppqz (P real and positive) carrying conformal weights {hpq,0) and (0, /i^^) (see Eq. ( ^.45| )), respectively. 
Moreover, {i'pqz,fpqz) transforms like a doublet (singlet) of the [/(l/l) Kac-Moody algebra in the holomorphic 
(antiholomorphic) sector, while {ippqzTfpqz) transforms like a doublet (singlet) of the [/(l/l) Kac-Moody algebra in 
the antiholomorphic (holomorphic) sector. Because (p_^) and {ip_,(p_) transform like doublets under diagonal 
[/(l/l), we must then have 

rV+(z,z)) ( i'pqziz)fq'zi^)] (^_{Z,Z)] ( ^pqz(z)fq,ziz)] , . 

U+(Z,Z)J ^ ['Ppqz{z)fq,,{z)j ' [^-{Z,Z)} ^ [ ^pqM f q' zi^) } ' ^ ' 

where ~ means obeying the same transformation law under diagonal [/(l/l). The real numbers p > 0,q, q' arc thus 
not a rbitrary, since the J^, j^^rj^^f]^ charges of the left-hand side must equal those of the right-hand side. From Eq. 
(3.56) and table |[ conservation of the current yields the constraint 



while conservation of the current yields 
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Current conservation fixes p = 1 and only leaves q = q' undetermined. We thus have found that 



f i^l, 1/'+, 'P+ ) 
[ tpi, ip^_, '4>_, J 



''I'lqzfqzi V'lqzfqzi i^lqzfqzi flqzfqz 
^ ^Iqzfqz^ ^Iqzfqz^ ^Iqzfqz^ ^Iqzfqz , 



(6.60) 



Again, ^ means obeying the same transformation law under diagonal C/(l/l). Having derived the correlation func- 
tions of the primary fields on the right-hand side of Eq. ( 6.6C ) , the correlation functions among {i/j^, , • • • , } can 
be calculated. But since {?/;^, • • • , are local operators, their correlation functions must be single-valued. The 
condition of locality then further limits the real valued q. For example, 



((p+(z,z)¥'U0'0)) 



is only single- valued if 



2h 



(6.61) 
(6.62) 



where I is an integer. In other words, using Eq. (6.45) we find that q is parameterized by an integer I 



q = Ik — 



1 



1 - 



I e Z. 



(6.63) 

, } and powers thereof are 



One can verify that Eq. (6.63) insures that all correlation functions among 
single- valued. 

Finally, we show that the constraint of locality on the OPE among ■ ■ ■ t^P-\ forces the Kac- Moody central 
charge k to be quantized. Without loss of generality, we only consider the OPE between ^^fz, z) and V3+(0, 0). Except 
for the residues, it is th e same as the OPE between ipiq^{z) and ^Piq^i^) due to Eq. ( |6.39 ). It must therefore be of 
the form (see Eq. ( S.51 )): 



(pAz,z)LpAQ,Q) 



(2g) 



.(0) 



(6.64) 



Locality again requires that 



'2(2?) 



2h 



fc' 



(6.65) 



be an integer, i.e., k must be the inverse of some integer m G Z. We thus see that the [/(!/!) critical point of C 
must be labelled by the triplet {I, m, A; ■) with m G Z and A; ■ being a real number. 



C/(l/l) 



E. Operator content of the t/(l/l) x (7(1/1) critical points 



We assume that we are at the critical point labelled by the triplet [l, m, k^). The full local operator content of this 
critical point can be obtained from the operator algebra generated by the currents (J, j, rj, r]) and the local operators 
{ip+j •••,'/'-}• In the following we want to calculate the conformal weights of some local composite operators. We 
begin with 



neZ. 



(6.66) 



With a calculation along the lines of the previous subsection, one verifies that the conformal weights {h , 0) of$, , 
are the same as the conformal weights (h^i^,, 0) of ^p,gi^, where p' = n,q' = nq. Hence, with the help of Eqs. ( |6.55 ) 
(|05|) and (|6J3|), 



h = —'m\n\ + —n'^{l — m), 7i G Z. 



(6.67) 



Similarly, has (0, h) for conformal weights. More generally, we introduce (compare with Eq. ( 3.30. 



V'+V'- , > 0,7l2 > 0, 

ni< 0,7^2 >0, 

V'+V-""'' 0,7^2 <0, 

^i^"V^"^ n,<0,n,<0, 



(6.68) 
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where ni,n2 are integers. The conformal weights {h,h) of <I>„ „ are equal to: (i) those of iPpiqi^ in the holomorphic 
sector with p' = \ni\,q' = (|ni| ± |n2|)(7, and (m) those o f ifi^i g'z antiholomorphic sector with p' = |n2|,(7' = 

{\n2\ ± \ni\)q. Thus, according to Eqs. (|^, ( |6^ ) and ( p^ ) 



11 = — m ni 
2 ' ^ 



1 



1"2 



1 

— rii nn 
2 ^ ^ 



I- 



4 



(6.69a) 
(6.69b) 



Hence, the scahng dimension A„j„2 ^ h + h and spin s,ii„2 = h — h of „ for Wj^, ^2 G Z are 



A„^„^ = ^m{\ni\ + In^l) + ^{nj + nl){l ~ m) + n^n^ 



(6.70a) 
(6.70b) 



respectively. Notice that we recover the conformal weights of Eq. (3.30b) if ~ 



We see from Eqs. ( |6. 67 , 3.70b) that the spins of and $niri2 ^-re either integers or half-integers depending 



on whether the labels Z,m S Z of the critical point are even or odd . Th e composite operators $+„,$_„, and 
are thus local for all integer valued n,ni and n2, according to Eq. (3.54). 

It should be noted that for any critical point (/, m, fcj), the condition that V'± and Lp± have half-integer spin (although 
not necessarily i) must hold. This is so because the transformation corresponding to a rotation in the (r, x) plane by 
27r changes the sign of the spinors ip^ and ip^, while it leaves the impurity potential unchanged. Hence, this property 
should be preserved at the interacting fixed point. It then follows from Eq. (6. 70b) that I is an odd integer. Other 
constraints on the possible critical points {l,m,kj) require additional symmetries. For example by assuming that the 



interaction C 



U{1/1) 



is invariant under rotations of the Euclidean plane (a;,r), it follows that only Z = 1 is allowed. 



Indeed, at the free theory critical point (1, 1,0), V'± a-nd ip± have spin ±i. By assumption, this quantum number is 



unchanged at the new critical point {l,m,kj). From Eq. (6.70b), the condition that, say sio = ^, thus requires I — 1. 

In summary, the possible critical points Cq + CfJ^^^^^ with U{1/1) x [/(l/l) symmetry are labelled by the odd 
integer I, the integer m and a real number kj: (/, tti, k^). Many correlation functions can be calculated exactly from 
the OPE. Although {l,m,kj) enumerates all possible critical points with currents obeying a [/(l/l) x [/(l/l) Kac- 
Moody algebra, it is not clear whether or not all these critical points can be realized by a theory with a Lagrangian 
solely constructed from the fcrmionic spinors tpi: j '0± ^^'^ bosonic spinors (p\. , Lp± . We were only able to identify 
the line {1,1, kj), kj > 0, with a Lagrangian realization of the random U{1) vector gauge-like model. Finally, the 
same approach can be used to identify all critical points with currents obeying the G x G Kac-Moody algebra, where 
G is a graded Lie algebra. This we have done explicitly for the compact group U{2/2) a nd th e non-compact group 
[/(1, 1/1, 1). In both cases, all critical points are labelled by the triplet {I, m, kj) and Eq. (6.69) holds. 



VII. DISCUSSION 



We have studied a model of massless Dirac fermions in two spatial dimensions and in the presence of static random 
vector potentials which are distributed according to a white noise probability distribution. We have shown that this 
model is an example of a non-unitary two-dimensional GET with a highly unusual operator content, since there exists 
an infinite hierarchy of negative dimensional operators. The physical origin of this hierarchy of negative dimensional 
operators is the fact that the single-particle Green function is log-normal distributed with respect to realizations of the 
randomness. In turn, the non-selfaveraging nature of the single-particle Green function is related to the multifractal 
nature of a typical wave function with vanishing energy, although it is still an open question of how to calculate the 
multifractal scaling exponents associated to the critical (i.e., vanishing energy) typical (i.e., for a fixed realization of 
the disorder) wave function from the field theory. We have argued that such disordered critical points are unstable 
in a very special way, if perturbations corresponding to, say, randomness in the mass are present in the lattice 
regularization of Dirac fermions. Indeed, in addition to the strength of the mass-like randomness, the shape of its 
probability distribution induces infinitely many relevant perturbations through all the cumulants of the probability 
distribution. 

On the one hand, it would appear from our study of disordered critical points with [/(l/l) x [/(l/l) graded 
Kac-Moody algebra that the existence of infinite hierarchies of negative dimensional operators is related to the 
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multifractality of a typical critical wave function, and leads to a very special instability of the critical point. On 
the other hand, there are numerical evidences that the plateau transition in the integer quantum Hall effect is also 
characterized by negative dimensional operators 0. If so, they clearly do not destabilize the disordered critical 
point corresponding to the plateau transition. The most obvious difference between the plateau transition and Dirac 
fermions with random vector potentials is that the averaged density of states is non-critical in the former. This is very 
suggestive of requiring spontaneous symmetry breaking of the chiral symmetry of disordered Dirac fermions (there is 
no Mermin- Wagner theorem preventing this from happening in the context of a non- unitary CFT). Hence, it might be 
too restrictive to assume as we did that there is no spontaneous symmetry breaking of the t/(l/l) graded symmetry. 

A disordered system related to Dirac fermions in the presence of random vector potentials is the two-dimensional 
XY model with quenched random bond coupling studied by Rubinstein et al. Indeed, one can use the equivalence 
between the massive Thirring model and the Sine-Gordon model | |4^ , on the one hand, and the equivalence between 
the Sine-Gordon model and the Coulomb gas pO| , on the other hand, to relate the random XY model to interacting 
Dirac fermions in the presence of random vector potentials. For a Gaussian distribution of the randomness with small 
enough strength (variance) < cr < o'er, Rubinstein et al. predicted a phase with algebraic decaying correlations at 
finite temperatures for the disordered system. Korshunov [l6l| has argued that this intermediate phase with quasi- 
long range order is not stable. The source of the instability in Korshunov's analysis is very reminiscent of the special 
instability induced by large powers of the mass operator for our Dirac spinors. It would thus be extremely useful to 
understand in the context of the replica trick how to describe the special instability of non-interacting Dirac fermions 
in the presence of vector gauge-like disorder. 
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APPENDIX A: U(l/1) SUGUWARA CONSTRUCTION 



a. Semi-direct algebra 



We assume the formal construction |^ 

T{z) = ^1 [J J + r^fi - nv) [z] + {J J) iz) (Al) 

for the holomorphic component of the energy momentum tensor out of the holomorphic components of the currents 
J, j, 77,77 obeying a [/(l/l) Kac-Moody algebra. Above definition is meaningless until what is meant by the products 
of currents at the same point is specified. To this end we define the Laurent series 



r(z) = ^z-"-i/,'j, r e{J,j,77,^}- 



(A2) 



The algebra obeyed by the normal modes of the currents I°-{z) is 

where the bracket [A, B]^ is the anticommutator {A^B} for A,B two fermionic operators, and the bracket [A, 

is the commutator [A, S] for all other cases. According to Eq. (6.29), the only non-vanishing central extensions k"-^ 

and structure constants /'^''^ are 



^2 ^ ^21 
p233 _ f323 



iLl _ i-zi _ 97 j^-zz _ 1 1,34 _ _J, 



22 _ u ;,34 _ ,„43 _ f,^ 

244 _ f424 _ I o -f341 _ /•431 



jzid _ __j3:i;3 _ 2 / = — / = +2 / = / = +1 



(A4) 



40 



From the algebraic point of view, k and can be taken to be arbitrary complex numbers. The 12 commutators and 
4 anticommutators for the generators of the C/(l/l) Kac- Moody algebra are explicitly given by 

[Jm^Jn] = 0, [JmJn] = 2fcm(5„_„, [J„,?7„] = 0, [J„,^„] = 0, (A5a) 
bm:4] = 2fcm(5„__„, [j,„,i„] = A:jm(5„_„, [i„,?7„] = -2?7„+„, [j™,f?„] = +277„+„, (A5b) 

['7m''^rJ— 0; [^m J Jn] = +2?7m+ru {'?m7^n} = 0i {^m i 'ynl = ~^"^'^m,-n + "^m+n ; (A5c) 
[V,n, Jn] = 0, [ri,nJn] = ''^Vrn+n: {Vm^Vn} = +km5,^ _,^ + J,„+„, {?7^,fy„} = 0. (A5d) 

Having defined the Kac- Moody algebra C/(l/l), we introduce the convolutions of the currents through 

-^jn ^1 ^ ^ • {."^m+m'j—m' Vm+m'V—m' ~ Vm+m'V—m') ■ '^^2 ^ ' ■ Jm+m' ^ —m' ■ 

m'G Z m'(^ Z 



1 {j-m'Jm+m' ~ V-ni'Vin+m' V-in'Vm+m') ^^2 J~m"^m+m'- {^^) 



m+m'>0 m+m'>0 

The parameters 2 complex and to be determined below. Here, the normal ordering is defined by the second 
line on the right-hand-side. We now show that the property 

[L^,a = -nC+„, Vm,ne Z, (A7) 

holds if and only if 

1 4 - A:, 

= 2fc' ""'^ -if- 



Notice that it follows from Eq. ( A7) that all are left invariant by the generators /q , a = 1, • • • , 4. 
PROOF: Let to, n be arbitrary integers. Decompose L„j into three contributions; 

(A9a) 

^,„+„' :, (A9b) 

(A9c) 



7-1 


= '^l H 


'^m-\-7n' j —m' ' ' 




m'e z 












m'e z 








rn-\-rn' '-^ —m' ' " 




m'e z 





The case a = 1 is straightforward. It yields a constraint on k.-^. Indeed, due to Eq. (A5a), [L^^ J^] = [-^^5 Jri[- Now, 



m+m'<0 



m+?n'>0 



= -n 2fcKi J„+„. (AlO) 
The case a = 2 is more involved. It yields a constraint on Ki and The first contribution is 

m+m' <0 

"I" '^l ^ ] {j-m'[Jm+m' t Jn] [j-m' ^ Jnl'^rn+m') 
m+m'>0 

= ^^111 ('^m+™'fcj(-"i')'^-m',-n + 2fc(m + to')(5,„+„, j„„.- ) 
m+m' <0 

+ Kl ^ (j_,^,2fc(TO-KTO')^m+m',-n + J^+m' ) 

m+m'>0 
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The second contribution is 



m+m'>0 



m+m' <0 

fm'<0 m+m'<0 

^ ^ iVn—m'Vm-\-m' ~^ Vn—m'Vm+m'^ ^ ^ (jl—m'^m-\-n-\-m' ^~ V—m'Vm-\-n-\-'m'^ 



-2ki 



m+m'>0 



m+m'>0 



(A12) 



It is very important to account for the fermionic sign in the normal ordering of L^. This allows us to reduce the 
infinite sum to a finite sum. To see this, introduce the integer 



I = min{— m, — m + n}, 
and change variables in the second sums in the square brackets: 

-2k 



m'<—m m" <—m-\-n 

\Vn—m"Vm+m" Vn—rn'Vm+m" ) 

m'> — m 'm"y — rn-\-n 

^ y (jlm+m' Vn—m' ~^ Vm+m'Vn—m'^ ~^ ^ ^ iVn—m'Vm+m' ~^ Vn—m'Vm+m', 



m'=I 

We are thus left with the finite sum 



m'=I 



/+|n|-l 



[-^m'in] ~ Inl ^ ({^m+m''^n-m'} + {'?m+m''^n-m'}) 

' ' m'=I 

' ' m'=I 

+ fc(m + m')(5^_,_^,__„_,_^, + Jm+n^ 



Finally, the third contribution is 



m+m'<0 



+ '^2 ('^-m'['^m+m'' in] + ['^-m''in]'^m+m') 

m+m'>0 

m+m'<0 

+ K2 E (•^-m'("^ + "j')'5m+m',-n + (-m')<5_„/,_„ 



m+m'>0 



(A13) 



(A14) 



(A15) 



(A16) 
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Collecting all three contributions yield 



(A17) 



Without loss of generality, we only consider the case a — 3, the case a — 4 being similar. We show that and 
K2 are not overdetermined. We only need the commutators of r/„ with and L^^^ since J,„ commutes with 77^. The 
first commutator is 



m-\-m' <0 



■m+m'>0 



m-i-rn' 



m+m' <0 



7n4-m'>0 



(AI8) 



m'e z 



The second commutator is 



(A19) 



(A20) 



m-\-m' <0 
m+m'>0 

+ 5Z (^m+m'^(^'™')'^-m',-n + '?m+m'^n-m' + '^("^ + "^')^m+m',-n^-m' + ^m+n+m' '7-m' ) 
m+77i' <0 

+«! (fc( — m )<5_m',-n^m+m' + Jn-m'Vm+m' + ^-m'^("^ + )'^)n+m',-n + V-m'J m+n) = 

m+m' >0 

m'e z 

Collecting all terms, we find 
This concludes the proof. 

b. Virasoro algebra 

We now show that the set {i,„},„g ^ generates a Virasoro algebra with vanishing central charge, provided Ki and 
are chosen as in Eq. ([A^): 

PROOF: Let m,n be arbitrary integers. The first commutator to be considered is 

n+n'<0 n+n'>0 

— '^ri+ri'Jm-n' ~ + )'^r?j+n+n'.?-n' 

n+n'<0 n+n'<0 



(A22) 



n+n'>0 



n+n'>0 
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With the help of the change of variables n' = m + n" in the first summation on the last two lines and with the 
introduction of 

/ = min{— n, — n — m}, (A23) 

we can rewrite the right-hand-side as 

[Lm,Ll] = Hi i^ + ^")Jn+m+n"j-n" ~ 1^1 + "^m+n+n' J -n' 

n+m+n" <0 77i-\-7i-\-7i' <0 

+ '^1 Y (™ + "")j-n"'^«+m+n" - «1 Y + -n' Jm+n+n' 



■n+m+n">0 m+n+n'>0 

I+\m\-l I+\m\-l 

I I ri'—I ' ' n'—I 

/+|-m|-l 

(m - n) i;^+„ - Ki — V -I- n')[J„,+„+„/, j_„/] 

TO 

' ' n' — I 

I+\m\-l 

m 



(to - n) i,i„+„ - 2fcKi — Y in + r>'')im + n + n')S,„^^^^, „^, 

I+\m\-l 



TO 

' ' n'=I 



= ("^-")im+« - — (5„^_„ ^ (n-Hn')n'. (A24) 

' ' n'=I 

Because the commutators of Jp with vanishes for any integers p and q, one immediately obtains from [L^, L^] that 

[L,„,L3] ^ ^3^^^^ (A25) 

The proof that the Virasoro central charge vanishes is then established by showing that 



[L^,Ll] - (TO-n)L2„+„ + 2fcKi ^<5„,_„ ^ (n -f (A26) 



n'=I 

To show this, start from 

[^rm^'n] ~ '^l [^mJ {Vn+n'V-n' ~ Vii+ii'V-n')] ^ '^1 [-^m: {V-7i'Vn+n' ~ V-n'Vii+i. 

n+n'<0 n+n'>0 



'^l {Vii+ii'Vm-ii' Vn+n'Vm-n') '^1 + ('?m+n+n'^-ra' Vm+n+n'V-7i') 

n+n' <0 n+n'<0 

-/ti ^ n' {nm-n'Vn+n' " '^m-n' ) + ^1 + i^^n'^m+n+n' ^ V-n' Vm+n+n') ■ (^27) 



n+n'>0 n+n'>0 



With the help of the substitution n' = m + n" in the first summation on the last two lines, we can again extract 
(to — The reminder is a finite sum over anticommutators due to the fermionic normal ordering: 



I+\m\-l 

TO 



[■^m:^n] = (™- »^)^m+n - '«1T— T Y (" + '^O ({Wn+n' ' '7-„' }-{??,„+„+„', 77-„' }) • (A28) 



TO, 



The term J^+n cancels from the difference of the two anticommutators and we are left with twice the contribution 
from the central charge k to the OPE between 77 and 77 as given in Eq. ( A26| ). This concludes the proof. 



APPENDIX B: LAGRANGIAN REALIZATION OF (7(1/1) x (7(1/1) CURRENT ALGEBRA 



In this appendix, we show that the ?7(1/1) x C/(l/l) Kac-Moody algebra of Eq. (A4) or Eq. ( |6.29 ) with real central 
charges {k,kj) is associative. To prove associativity, it is sufficient to construct a Lagrangian realization of a free 



44 



field theory, whose currents obey the algebra of Eq. ( 6.29| ) with {k,kj) = (1, fcj), kj > 0. Associativity for arbitrary 
values of {k,kj), where k 0, then follows on purely algebraic grounds by independent rescalings of the currents 
(J, j, rj, ij). We have already seen a Lagrangian realization of the current algebra with central charges {k, kj) — (1,0) 
in Eq. (6.1). We now construct a Lagrangian realization of [/(l/l) x U{1/1) Kac-Moody algebra with central charges 
{k, kj) = (1, fcj), kj > 0. (For a similar trick, see Ref. |6^]) Consider the partition function 



A.. 



(dj) (9,0) . 



(Bla) 
(Bib) 
(Blc) 



Defin e the holomorphic and antiholomorphic components of the free real scalar field (j){z, z) by (compare with Eq. 

(Hi)) 



(z,z) = </.(z) + 0(z), (0(^)0(0)) =-2 Inz, (^(z)^(O)) --21nz, (0(z)^(O)) = 0. 



(B2) 



The real scalar field ((> is taken to commute with all spinors. Without loss of generality, we only consider the 
holomorphic sector. Define the holomorphic component of the currents 



Jiz) = {4,l^P+) (z) + (^;^+) (z), 

J-(^) = {z) 
r]{z) = (ifi'ip+'j (z), 



f d^iz), 



(B3a) 

(B3b) 
(B3c) 
(B3d) 



The structure constants {/^fcj,} ^'^^ ^^'^ OPE of these currents are those of J7(l/1), since cj) commutes with the spinors. 
By assuming that currents do not acquire vacuum expectation values, the only change in the algebra takes place in 
the OPE j{z)j{Q) = Note the similarity between Eq. (Bl) and the effective critical theory defined by Scr' in 
Eq. (EH). 
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